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We review and compare theoretically and phenomenologically a number of possible 
family symmetries, which when combined with unification, could be important in 
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non symmetric Yukawa matrices. Grand Unification or not. Our main focus is on 
anomaly-free U{1) family symmetry combined with SU{b) unification, although we 
also discuss other possibilities. We provide a detailed phenomenological fit of the 
fermion masses and mixings for several examples, and discuss the supersymmetric 
flavour issues in such theories, including a detailed analysis of lepton flavour violation. 
We show that it is not possible to quantitatively and decisively discriminate between 
these different theoretical possibilities at the present time. 
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1 Introduction 



The hierarchy of quark and charged lepton masses and the small quark mixing angles 
has been one of the most puzzling aspects left unresolved by the Standard Model. 
The recent discovery of neutrino masses and mixings has provided further clues in the 
search for the new physics Beyond the Standard Model which must be responsible 
for the pattern of fermion masses and mixing angles. One promising approach to 
understanding the fermion spectrum is the idea of family symmetry, and in particular 
the idea of a family symmetry as originally proposed by Proggatt and Nielsen [1] . 
Such an approach was given considerable impetus by the observation that in many 
string constructions additional U{1) symmetries are ubiquitous, and furthermore such 
a gauged broken U{1) could provide a phenomenologically viable candidate family 
symmetry by virtue of the Green- Schwartz anomaly cancellation mechanism [3] which 
provides a string solution to the no-go theorem that anomaly freedom requires such 
symmetries to be family independent [2] . As a result of this a considerable literature 
has developed in recent years based on string- inspired U{1) family symmetries [4,5]. 

Many non-abelian family symmetries have also been considered, for example based 
on SU{3) family symmetry [6], and also textures and analyses of fermion masses 
have been done not using any family symmetry. At the present time some very 
successful approaches exist, and others that may with modification also be effective. 
Family symmetries can be abelian or non-abelian, they can require symmetric Yukawa 
matrices or not, they can be imposed with or without an associated grand unified 
theory, and so on. Criteria that could be used to choose among possible approaches 
include not only describing the quark masses and mixings, and the charged lepton 
masses, but also neutrino masses and mixings, supersymmetry soft breaking effects 
(since particularly the trilinear couplings are affected by the Yukawa couplings) , how 
many parameters are used to describe the data, whether some results such as the 
Cabibbo angle are generic or fitted, and more. One of our main goals here is to 
look at the various possibilities systematically and see if some seem to be favoured 
by how well they do on a set of criteria such as the above listed ones. Presumably 
family symmetries originate in string theories, and are different for different string 
constructions that lead to a description of nature, so identifying a unique family 
symmetry (or a subset of possible ones) could point strongly toward a class of string 
theories and away from other classes. At the present time this approach is not very 
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powerful, though it gives some interesting insights, but better analyses and additional 
data may improve it. 

In this paper we shall consider U{1) family symmetries and unification as a viable 
framework for quark and lepton masses and mixing angles in the light of neutrino 
mass and mixing data [7], using sequential right-hand neutrino dominance [8] as a 
guide to constructing hierarchical neutrino mass models with bi-large mixing. As 
has been pointed earlier [9], models which satisfy the Gatto-Sartori-Tonin relations 
(GST [10])^ require the presence of both positive and negative Abelian charges. As 
we will discuss, the sequential dominance conditions require also the presence of both 
positive and negative Abelian charges, and hence at least two flavon fields of equal 
and opposite charges. These models however result in complicated U{1) charges, on 
the other hand Non-GST models have a simpler charge structure and may be possible 
to realize in a more general context. In this work we also consider non GST cases. 

We shall consider U{1) family symmetry combined with unified gauge groups 
based on SU{5) and 5*0(10), assuming a Georgi-Jarlskog relation, and also consider 
non-unified models without such a relation. We will present new classes of solutions 
to the anomaly cancellation conditions and perform phenomenological fits, and we 
will compare the different classes of U{1) to each other and to non- Abelian fam- 
ily symmetry models based on SU{3) [6], by performing specific phenomenological 
fits to the undetermined coefficients of the operators. Finally we will consider the 
implications of such an approach on flavour- changing processes in the framework of 
super symmetry, leaving a detailed analysis for a future reference. 

The layout of the paper is as follows. In Section |21 we consider the general condi- 
tions for Green-Schwartz anomaly cancellation, and move on to describe the classes 
of solutions, by whether they are consistent with SU{5), 5'O(10), Pati-Salam unifica- 
tion of representations, generalized non-unified relations, or not at all consistent with 
unification. Having found these solutions, we move on in section El to re-parametrize 
in terms of differences in U{1)f charges. In section 0] we consider the constraints on 
the Yukawa textures from requiring acceptable quark mixings and quark and lepton 
masses. Then in section El the constraints from getting acceptable neutrino masses 
and mixings from single right-handed neutrino dominance (SRHND) models, which 
are a class of see-saw models. In section we construct solutions which are consis- 
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tent with SU{5) unification, the Gatto-Satori-Tonin (GST) relation [10], and correct 
fermion masses and mixings. In section [7| we construct solutions which are consis- 
tent with SU{5) unification, correct fermion masses and mixing angles but which are 
not consistent with the GST relation. In section |H1 we construct solutions which are 
not consistent with SU (5) unification. In section El we take some of the solutions 
constructed in section El and section [71 and fit the arbitrary 0(1) parameters to try 
to closely predict the observed fermion masses and mixing angles. Then in section 
ITUl we briefiy consider whether fiavour changing processes will be dangerously high 
in these models, presenting two specific scenarios: a non minimal sugra possibility 
and a string-inspired mSUGRA-like scenario which is expected to be (or be close to) 
the best-case scenario for fiavour-changing and for which we check explicitly — 67 
Finally, we conclude in section ITTl 

2 Anomaly Constraints on U(l) Family symmetries 
2.1 Green-Schwartz anomaly cancellation 

Consider an arbitrary U{1) symmetry which extends the Standard Model gauge group. 
If we were to insist that it does not contribute to mixed anomalies with the Standard 
Model, we would find that the generators of U{1) would be a linear combination of 
Weak hypercharge and B — L [2]. This clearly is not useful for family symmetries, so 
we need to use a more sophisticated way of removing the anomalies, Green-Schwartz 
anomaly cancellation [3]. In this case, we can cancel the mixed U (1) — SU (3) — SU [3) , 
U{1) - SU{2) - SU{2) and f/(l) - f/(l)y - f/(l)y anomalies, A3, A2, and Ai if they 
appear in the ratio: 



where we have included the relations to the anomalies of the anomalous fiavour groups 
and the gravitational anomaly; ki are the Kac- Moody levels of the gauge groups, 
defined by the GUT-scale relation: 



If we work with a GUT that has the canonical GUT normalization, we find: 



A3 : : : Au(i) : Ac = h : k2 : h : 3kuii) : 24, 



(1) 



glh = glk2 = glki 
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But we still require that the f/(l) — f/(l) — U{1)y anomaly, A[ vanishes. Now, the 
anomalies are given by: 

A, = iTr[{r«,T«}T;]. (4) 

We then use the fact that {T„, = 6abl for SU{N) and {Y, Y} = 2Y^ for U{1)y to 
obtain: 



A, 



A, 



^{2qi + Ui + di 

1=1 

3 

^(3gi + k) +K + K 
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E.Oj Suj 2 , 3/j 6ej , 3 



(5) 
(6) 
(7) 



= Y.^-i + ^<-dl + lj-el) + {hl-hl) = Q 



i=l 



Since in the mixed anomalies of the U{1) group with the SM gauge group that cancel 



Field 


Qi 




Li Ei 


Ni 


Hu 


Ha 


Charge 








rii 


K 


hd 



Table 1: Fields and family charges 

via the Green-Schwartz mechanism wherever a charge appears, it appears in a sum, 
we parameterize the sums as follows [12]: 



X + u, ^^Mj = X + 2u, 
1=1 

3 

^di = y + v, 

i=l 
3 



1=1 
3 



i=l 



(9) 
(10) 



i=l 



X, (11) 

-z, hd = z + {u + v). (12) 
Substituting Eq. ©-Eq. ^ into Eq. ®-Eq. © we find that they satisfy Eq. (0): 

(13) 
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3 1 

Az = A2 = -Ai = -[?,x + Au + y + v\, 



which shows that the parameterization is consistent. However we need to find those 
solutions which also satisfy A[ = 0. We will see how we can achieve this for different 
cases. Since the proposal of the GS anomaly mechanism it has been known that 
the easiest solution, u = v = 0, leads to a SU (5) or Pati-Salam group realization 
of mass matrices. Another possible solution is to have u = —v ^ 0. Both these 
forms admit a SUSY /i term in the tree level superpotential at the gravitational 
scale. However given the form of Eq. (jHll-Eq. (|T2|) one can try to use the flavour 
symmetry in order to forbid this term, allowing it just in the Kahler potential and 
thus invoking the Giudice-Masiero [11] mechanism in order to generate the /i of the 
desired phenomenological order. Therefore apart from the cases w + f = we examine 
plausible cases for u ^ —v ^ ^. Of course in the cases -u = f = 0,^ = —v ^ one 
can use another symmetry to forbid the /i term in the superpotential, however it is 
appealing if the flavour symmetry forbids the ^ term at high scales. 



2.2 Anomaly free A'-^ with u = v = solutions 

In this case the parameterization simplifles and in fact we can decompose the U{1) 
charges in flavour independent and flavour dependent parts 

u = \f+n- (14) 

The flrst term is flavour independent because it just depends on the total sum of the 
individual charges and the // are flavour dependent charges. We can always flnd x 
and y which satisfy 

3 

1=1 

In this way A'-^ can be expressed in flavour independent plus flavour dependent terms 

^'i = Afi + Afd- (16) 

Following this, with the unfortunate notation that we have a new tt, completely 
unrelated to the u that we have already set to zero, we then have: 

A!^ = A[pj + A[pj^ 

= \[-q' + 2u'-d' + l^-e']+j2i-<ll' + H'-<' + l'^-e',') (17) 

^=l 
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Now it is clear that the terms in the square bracket in Eq. (fTTj) are family independent. 
It turns out that the square bracket term is automatically zero in this case, since from 
EqsinilTH we have: q = u = e = x and I = d = y. Then we have to make the family 
dependent part (the second term in Eq. (fTTjl ) vanish. 



2.2.1 5f/(5) and 50(10) type cases 



One way to make the family dependent part vanish, A'^p^, = , is to set /j = di and 
qi = Ui = Ci ^. This condition would be automatic in S'f/(5), but in general such 
a condition on the charges does not necessarily imply a field theory SU (5) GUT to 
actually be present, although it may be. 

Since the generic Yukawa structure is of the form: 



gl/i+gi+'i/l 
^\fi+q'i+hf\ 



^\f2+qi+hf\ 

gl/2+l22+/l/ 
Af2+q-i+hf\ 



gl/3+gi+/i/l 

gl/3+'?2+/l/| 

Ah+qa+hf\ 



it is clear that the SU{J)) relations di = k, qi 
the form: 



Ui 



Ci lead to Yukawa textures of 



y-b 



Y" 



Y' 



|2ei-2e3| ^\e-i+e2-2e-j,\ ^\ei-e-j,\ 



g|ei+e2-2e3 
glei-esi 



g|2e2-2e3| g|e2-e3[ 



r|e2-e3| 



.|0| 



\ll+ei+ha\ g|«2+ei+/id| gl's+ei+^d 
\h+e2+hd\ f:\h+e2+hd\ ,|i3+e2+^d 



Y 



e 

d T 



h+es+h^l Ah+es+hdl Jh+es+hdl 



(19) 



(20) 



(21) 



Note that the up matrix is approximately symmetric, due to the assumed SU{5) 
relation of charges. The reason why the textures above are approximate is that 
each entry in each matrix contains an undetermined order unity flavour dependent 
coefficient, generically denoted as a{j = 0(1). We shall continue to suppress such 
coefficients in order to make the discussion less cumbersome, but will return to this 
question when we discuss the numerical fits later in the paper. We have also assumed 

^The reason that the charges are unprimed here is that if it is true for the primed charges, it is 
also true for the unprimed charges 
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that the up and down Yukawa matrices are described by a single expansion parameter 
e. The possibihty of having two different expansion parameters, one for the up sector 
and one for the down sector, will also be discussed later in the paper. In order to 
have an acceptable top quark mass, we have required that /i„ + 2e3 = 0, in which case 
the smallness of the bottom quark mass can be due to + + I3 ^ 0, and we are 
free to have a small tan /3, because we don't need large tan /3 to explain the ratio ^ 
on its own. 

Also note that, as expected from the SU{5) relation of charges, the down and 
electron textures are the approximate transposes of each other, ^ O^'^)'^- Such 
a relation implies bad mass relations for between the down type quarks and charged 
leptons, but may be remedied by using Clebsch factors such as a Georgi-Jarlskog 
factor of 3 in the (2,2) position of the charged lepton Yukawa matrix. 

If we were to look at the case x = y, then we would have a solution suggestive of 
unified 5*0(10) GUT symmetry, for which U = qi = Ui = di = Ci. The same comments 
above also apply here, namely that such a condition on the charges, though consistent 
with an 5*0 (10) GUT does not necessarily imply a field theory realization of it. The 
matrices Eq. (jl9|) -Eq. (j20|) would all become equal to the same symmetric texture in 
EqUni in the S'O(IO) case that x = y. 



2.2.2 Pati-Salam type cases 



In this case, applying the Pati-Salam constraints on the charges, 

k = , Ui = di = e 



(22) 



so we can immediately see that also for this choice of charges both the the flavour 
independent and dependent parts in Eq. (fT7|) vanishes. We have also included the 
right-handed neutrino charges, which do not enter into the anomaly cancellation 
conditions, Eq. ©-Eq. (jH)), but with a Pati-Salam group should obey the relation of 
Eq. ()22|) . Thus in this case all the mass matrices have the form 

^ g|/i+ei+/i/| ^\h+e2+hf\ ^{h+e^+hfl \ 



Yf 



(23) 



y ^\l3+ei+hf\ ^\l3+e2+hf\ 

for hf = hu, hd- In this case we always need to satisfy x = ?/, in contrast with the 
generic case of S'f/(5) where it is not necessary x = y. So we can put one of the 
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charges in terms of the other two and the parameters x = y 



ei = x - (ea + 63), li=x-{l2 + h), ^ ei + 62 + = h + k + h- (24) 

We have already noted that the Pati-Salam constraints on the charges imply that the 
anomaly A[ automatically vanishes. It is also a remarkable fact that the constraints 
in Eq|21]do not in practice lead to any physical constraints on the form of the Yukawa 
texture in Eg 1231 In practice, assuming only that u + v = 0, one can start with any set 
of charges k, Cj which lead to any desired Yukawa texture, where the charges do not 
satisfy the anomaly free constraint in Eql24l Then from any set of non-anomaly-free 
charges one can construct a set of anomaly-free charges which do satisfy Eq l24[ but 
do not change the form of the Yukawa matrix in Eq l^ by simply making an equal 
and opposite flavour-independent shift on the charges as follows [30]: Cj — > Cj + A, 
h ^ h ^ In this paper we shall not consider the Pati-Salam approach in detail. 

2.3 Solutions with anomaly free A[ with u + v = {u,v ^ 0) 

In this case, we can repeat the analysis of the previous subsection, but with the 
general constraints. Note however, that since u + v = 0, = —z and hd = +z. 

Then we are left with the result that 



Having done this, we may substitute Eq. (j26|) into Eqs. (P- fT2|) Then we find that: 



A[ = ^ [6u' + 6xu + 2yu] - " '^<' + ' + e[') ■ 



(25) 



Note that the family independent part will vanish if 




(26) 




i=l 




(27) 



i=l 



i=l 



3 




(28) 



1=1 
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2.3.1 Yukawa textures for a sample solution 



At this point, we note tliat tliere will be a large number of solutions. However, one 
class of solutions that will easily be satisfied will be: 



,21 



Ui 



4L 



The same equation will hold for the primed charges: 



I' 



,21' 



Al'- 



(29) 



(30) 



We can now put Eq. (j3U|) into the anomaly, Eq. (|25|). In this case we find that: 

2 



x\Q-Q) + -y\l-l)+xy{A-2-2) 



i=i ^ 



8- 16 + 9) + e;'(2- 1 - 1 



0. 



(31) 



So we see that for this particular relation of leptonic and quark charges, we are 
automatically anomaly-free. 

Again, we see that, just as for the ti = f = case, we can specify everything 
by the leptonic charges U and Cj. However, in this case we will get three different 
textures. Specifically, we will get: 



ye 



rKi+ei+?i„| 
i(/i+2/2)+e2+/ii, 



g|i(«2+2ii)+ei+/i„| g||(Z3+2«i)+ei+/i.„I 



i(/i+2/2)+e3+?»i, 



U 



^\^il3+2l2)+e2+hn 
Ah+e3+hu\ 



.|/i+ei— hu 



Ji(-«i+4i2)+e2-/iu| g|§(-«l+4i3)+e3-/iu| 



g||(-«i+4«3)+e3-h„| g||(-«2+4i3)+e3-/i„| 



^\l(-l2+'il-i)+e2-hu\ 
^\h+e3-hu\ 



g|/2+ei-/iii 
c:\h+ei—hu\ 



\h+e2—hu 
\l2+e2-hu\ 
\l3+e2-hu\ 



\l2+e3~hu\ 
\l3+e3~hu\ 



(32) 



(33) 



(34) 



We note that this is a rather predictive scheme; we require that the diagonal elements 
are of the same order in the between the down and electron Yukawa matrices con- 
strained by the anomalies. Also, we require (at the very least) /s + 63 + /i„ = to get 
a correct top quark mass. 
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2.4 Anomaly free A[ with u + v ^ solutions 

In this case we can not decompose the expression of A[ into flavour independent 
and flavour dependent parts, but we can use for example the relation fif = 
Ef! + 2(/i(/2 + /s) + /2/3) such that we have 

A[ = -2{4u' + u{v + 3x + z)+ v{z - y)) -2^ (7^(/i(/2 + h) + /s/s), (35) 

f=u,d,l,e,q 

where gj = 1,-2,1,-1,1 respectively for / = q,u,d,l,e. However it is difficult to 
depart from here in order to find some ansatz which cancels the A[ anomaly. Instead 
we can generalize the kind of relations which in the limit of m = f = would give the 
SU{5) cases or the Pati-Salam cases. 

2.4.1 An extended SU{5) case 

Here a non-GUT case is considered, taken by generalizing the SU (5) relation between 
the charges. In the SU{5) case, we had qi = Ui = Cj and di = li. If instead we have 
the linear relations: 

gj = + a = Cj + 7, di = li + P, (36) 

From the parameterization of Eqs. (jSHHl), we see that in the limit of the u = v = 
we recover the SU (5) case. In agreement with the cancellation of anomalies then one 
should have 

Qi = Ui - - = ei + -, di = li + -. (37) 

In the expression of the A[ anomaly, as given in Eq. (jHl), the sums of squared charges 
cancel and we can write it just in terms of sum of charges, which we have parameter- 
ized in terms of u, v, x, y, 

V? 2 V 
a; = -10— - + 2u{x + v) + 2y- - 2z{u + v) = Q. (38) 

o o o 
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Thus we need to satisfy this equation in order to have anomaly free solutions. Re- 
quiring the condition of 0(1) top coupling we have 



K 




-2e3 - M, 


hd 


= 2u + v 


+ 263, 




1 ^2 ~t~ 


-2e3|, 




1 ^2 ~l~ 


7u Av , 


+ 263 H \ 

^3 3 ' 




= Ih + Cj 


+ 263 + 2M + t;|, 



(39) 

where C{Ylj) denotes the power of e for the {i,j) element of the correspondent Yukawa 
matrix. Note that although we did not begin with an a priori condition of having 
symmetric, the requirement of the 0(1) top coupling cancels the parameter u in all 
the entries of and so we end up with a symmetric matrix. 

2.4.2 An extended Pati-Salam case 

Following the extended SU{5) case, we look for solutions which in the u = v = 
limit reproduce the Pati-Salam case, so we should have the relations 

qi = k + a, Ui = di + (3. (40) 

Also 6j and rij need to be related to Ui by a constant, as in Eq. ()40|) . In these case in 

order to satisfy the G-S anomaly conditions we need 

u+{x-y) 2u v + {y-x) 
qi = , Ui = ei + —, di = Ci^ . (41) 

Thus the expression for the A[ anomaly is 
2 

A[ = --[Su"^ + Av^ + u{9v+llx-2y) + 2{x-yy -v{2x + y)] 

-2z{u + v), (42) 

and finally requiring the condition of 0(1) top Yukawa coupling we have 

7 /r ~ y \ 

hu = -z = -{k + e3 + u + 



hd = ^3 + ^3 + 2m + f + 
C(F,^) = |/.-/3 + e,-63|. 



3 

x-y 



^,^.ri^ , 4t> + 7m + (x — v) Av , 
C{Y^^) = \U + e,+h + e, + ^ ^ + yl' 

CiY^'j) = \li + e, + 2e3 + 2u + v\. (43) 
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3 A useful phenomenological parameterization 



So far we have discussed the anomaly cancellation conditions in f/(l) family symmetry 
models, and some of the possible solutions to these conditions, including some new 
solutions not previously discussed in the literature. It turns out however that the 
anomaly free charges themselves do not provide the most convenient parameters for 
discussing the phenomenological constraints on the Yukawa matrices arising from the 
quark and lepton spectrum. It is more convenient to introduce a new parameterization 
for the Yukawa matrices as follows: 



Yf 



(44) 



where / = u^d^e^v^ and we have introduced the parameters rf,r'jr,Sf,s'jr,kf which 
are defined in terms of the charges in Table 1 as: 

rj = f2- h r'j = fi- /a = q3 + U3 + K 
Su4 = 92-^3 = qi-q3 kd = qs + ds + hd 
Se,v = k- h = h-h ke = I3 + 63 + hd 

ku = h + n3 + K (45) 

In order to get an acceptable top quark mass, we require that k^ = 0. Note that the 
parametrization above is completely general, there is no information loss from the 
form of Eq. ()18|1 . and thus far we have not imposed any constraints on the charges 
arising from either anomaly cancellation or from GUTs. We now consider the sim- 
plifications which arise in the new parametrization when the charges are constrained 
by considerations of anomaly cancellation and GUTs, as discussed in the previous 
section. 



Simplification in SU{5) type case 

Consider the case where the family charges are consistent with the representations in 
an SU{5) GUT, di = k, and qi = Ui = e,: 

ke = kd Su,d = 'l"u,e S^^^ = g 

Se,u = Td s'^^^ = r'^ (46) 
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In this case, all of the parameters can be expressed purely in terms of the lepton 
charges: 

Su,d = ru,e = 62-63 s'^a = <,e = 61-63 

Se,u = rd = l2-h s'^,v = r'd = k- h (47) 

Note that this leads directly to the fact that (y^)^. The equality is broken 

by the arbitrary 0(1) coefficients. As discussed, the SU{b) charge conditions are 
sufficient to guarantee anomaly cancellation for the case u = v = 

Simplification in the extended S'f/(5) case 

In the case u + v 7^ 0, anomalies can again be cancelled by assuming the charge 
conditions in Eq. ()36p. If we take Eq. (jHBj) . we can again simplify Eq. In this 

case we find: 



Se,u = Td s'^^y = r'a (48) 

In this case we have that the texture of can be attained from by replacing 
kd with ke and then transposing. 

Simplification in the Pati-Salam case 

In the case of having charge relations consistent with a Pati-Salam theory, qi = U and 
Ui = di = Ci = rii, we can simplify: 

ke kd Su d ^e,u d i/ 

ku = ku ru = rd = re = r'^ = r'^ = r'^ = (49) 



4 Quark masses and mixings in SU{5) 

In this section we shall provide some constraints on the phenomenological parameters 
introduced in the last section, arising from the quark masses and mixings, assuming 
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the simplification in the SU (5) type case mentioned above. In SU (5) Eqs. 
imply the quark Yukawa matrices are explicitly of the form: 



e 



\2s' 



e 



\s'+s\ 



where we have written s 



e\s\ 



e\s\ 



yd 



e 



£:|s'+fc<il \ 

£-|s + fc£j| 



(50) 



Note that we are assuming a 



single expansion parameter e, and are suppressing 0(1) coefficients. Clebsch factors 
are also not considered, and only leading order operators are discussed. 

In order to determine the possible solutions for s, s', r^, and which success- 
fully reproduce quark masses and mixings one can numerically diagonalize Yukawa 
matrices and obtain the CKM matrix. However, in order to understand the behaivour 
of this structure it is quite useful to use the technique of diagonalization by blocks in 
the (2, 3), (1, 3) and (1, 2) sectors ^. The results are presented in the next subsections. 



4.1 Quark Masses 

Barring accidental cancellations the down quark Yukawa matrix may be diago- 
nalized, leading to the following eigenvalues: 

- 4(d^l + a22£2|r+fc|-|fcle-2^(/5|'-'3f)) ^ 

x(ai3£l«'+'=l+ a23ai2£l^'+'l+l^'+'^+'l-l'le-2^(^"-'^^')) + 
-(ai2£l"'+"+'^l-a32ai3£l"+'=l+l"'+^-|-|'^l)(a2i^rl"+"'+'^l-a23a3i£l'+''l+l''+^l"l^'l) 
(a22£l"+^+'=l - a2zaz2e\'^^\^\'-^^\-\^'\)e-^^f^'^-^^^ ' 

2/2- C2^3 (a22£l^+'^+'l - a23a32£l^+'l+l^+'|-|'l) -=^^(^"-^2'') 



e 



where we have suppressed the index d in order to make clearer the notation and re- 
scaled all the (complex) coefficients by 1/033, so that instead of having 033 we have 

■^Note that the extended SUih) anomaly free solutions examined in section 1^.4.11 leave the pa- 
rameters s, s', r^, kd invariant, as is clear by comparing Easl46land 1481 Hence the results in this 
section for the quark sector apply not only to the SU{5) type case but also the extended SU{5) 
anomaly free cases. 

^This only works if there is an appropriate hierarchy among the elements 
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1. Note that the down quark masses are given by: mf = yfvd/V^. Analogous results 
also apply to the up quark sector, with the replacements r ^ s, r' ^ s' , k ^ 0. The 
phases /3f correspond to the diagonalization matrices of the Yukawa matrices, whose 
notation is given in Appendix 

It is important to remark that in the case of positive charges all the elements of 
the first row of the Yukawa matrix contribute at the same order, s' + r' + k, to their 
correspondent lightest eigenvalue, so in these cases it is not possible to have the Gatto- 
Sartori-Tonin (GST) relation. However in the cases of having s and s' (analogous for 
r and r') with different sign, as in the example of Eq. (j58|) . we can have a cancellation 
in powers of e to the contribution to yi coming from the diagonalization in the (1,2) 
sector, which is the third term in the expression for yi in Eq. (jS^- On the other hand 
we can have an enhancement in the power of e of the contributions from the (1, 1) 
entry and the rotation in the (1,3) sectors, which correspond to the first and second 
term of yi, respectively, in Eq. This together with the condition C(F2i) = CiXu) 
are the requirements to achieve the GST relation. We will present examples satisfying 
and not satisfying the GST relation. 

We remark here the constraints from the bottom mass are 



since = 0((if„)) and tan/5 = {Hu) / {H^ . Thus in terms of charges we have 
hu = -{qs + Us) and hd = qz + Ms, for m = = 0, = 2g3 + 4 + u^. 

4.2 Quark Mixings 

We can also obtain the mixing angles in this approximation and compare to the 
required experimental values (see Appendix El- The mixing angles in the down 
sector, again dropping flavour indices, follows: 



rrib tan /5 = e' '''mj. 



kd = qs + ds + hd 



(52) 




t 



R 
23 



e*(/32^-/3f)a23£l^+'=l-l^l 



+ a23a22£^ 



\s+r+k\ + \s+k\-2\k\^iS,l^ 
.\s+r+k\ + \s+k\-2\k\^iili 
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n c-k'+^l 1/7 As+k\ + \s+r'+k\-\k\^2i(P^-f3{) , 

fR = ^3lg' '+a23Q2lg e ^2 \ A l U2 |^2|r+fc|-2|fc| 



''12 



R 

12 

6 



{a22e\'+'-+''\ - a23a32el^+'=l+l''+'=l-l'=l) 

(a22£l*+'"+'=l - a23a32el*+'=l+l'^+'=l-l'=l) 
-{^2 - Pi) - '^{^2 - /3f ), in = -{^^ - /?f ) - 2(/52^ - /?f ). (53) 

Analogous results also apply to the up quark sector, with the replacements — > s, 
— > s', fcrf — > 0. Note that in the case of positive s, s', r, r' and fc, the angles tfg ^'^'^ 
^23, of the left sector do not depend on r^, r'^, so they are equal, at first approximation, 
for the up and down sectors. Having the tangent of the angles expressed in terms of 
the Yukawa elements we can see directly their contributions to the CKM elements 
(VcKM = L'^L^^ in the notation of Appendix 

IT/ J c;"? _ pi{^l-'l>2)\ 

\^ub\ I"'12"'23 ■'13'^ I 



|4l 

\vu\ kf24-4e'(*^^l 



0.09 ~ (A^A) 



A 



\v,s\ |4| 

= |sf2-s^2e**'| = A^O.224 
Im{J} = 4(44s5'2sin($i)-4(«f2sm($2))-sr2sin($2-$i)), (54) 

with sfj = \sfj — e'^^^^s'^j\. The phases $i, $2 and ^Xij depending on the contributions 
that the mixing angles receive from the different elements of the Yukawa matrix and 
have a different expression in terms of the phases of the Yukaw matrix for different 
cases. For example when the elements (1, 2) and (1, 3) are of the same order and the 
right handed mixing angle in the (2, 3) sector is large, the $2 phase will be 



$2 = Arg 



'Vd I \rd4.R 
^12 ' ■'^ 13''23 

-'33 "T -'23''23 



(55) 



As we can see from the expressions in Eq. fl54j) involving $1, this can be associated 
to the U sector. When all the signalization angles in this sector are small, then this 
phase takes the form 

$1 = 012 - 022 ( 56) 
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where 0i2 and 022 are the phases of the 1^2 ^22 elements. Finally the phases 
which appear in s^, can be associated either with the U or with the D sector. 



U{1) relations 



Constraint 


Reason 


U{1) relations Constraint 


Reason 




"'23 


^|3<?3+<i3| ~(1,A3) 


nib 


~ A 


■^13 
^12 


^\s+rd+kd\-\ka\ ^ ^ _)^3^ 
^\s'+r'a+ka\-\kd\ ^ (_X4^ _)^5) 
^|2s+fcd|-|fcdl ^ _\4 
^|2s'+fcd|-|fcd| y ^6 


rns 
rrib 
"Id 
nis 
rric 
mt 
rriu 
rric 



^\s+kd\-\kd\ 
g\s'+kd\-\kd\ 
^\s'+rd+kd\-\s+rd+kd\ 



Table 2: Constraints on the parameters s, s' , r^, r'^ and from quark mixing angles and 
mass ratios. For the mixing angles we need to satisfy the conditions for up or down sector, 
where the analogous conditions for the up sector are obtained by making the replacements 
rd s, r'^ ^ s', kd — > 0. They do not need to be satisfied for both as long as for the sector 
in which they are not satisfied they do not give a bigger contribution than the indicated 
power. 



With the requirements of Table Q and the values of quark masses in Appendix 
, we can identify the viable solutions in the quark sector. One solution which has 
been widely explored is the up-down symmetric case for which we have x = y thus, 
fi = li = Ui = Ci = di = li. In this case hu = — 2e3 = —hd so ku = 0, kd = h = Ae^, 
but in this case we need two expansion parameters Eu and Sd to reproduce appropriate 
mass ratios and mixings, thus we have 



Yf 



|2s'+fe/| \s+s'+kf\ 
\s+s'+kf\ \2s+kf\ 

^\s' + kf\ ^\s + kf\ 



Js' + kf\ 

■f 

_\s + kf\ 



(57) 



We can think of fixing s + s', and then check for which choice of s we have appropriate 
phenomenological solutions. For example if we take s + s' = ±3 and 63 = {kf = 0, 
V/) we have 



1 6-2/2 1 



e 



|3| 



/ 

J3-f2\ 

L ^/ 



1 2/2 1 



^|3-/2| 

^1/21 
1 



The viable phenomenological fit for the case of quarks is for /2 = — 1 and /i 
/2 = 1 and /i = —4 [17]. In this case we have then x = y = ±3 respectively. 



(58) 



4 or 



17 



5 Neutrino masses and mixings in SRHND 



In this section we apply the requirements of getting acceptable neutrino masses and 
mixings by using a class seesaw model where h = h- These are a subset of a class of 
seesaw models called single right-handed neutrino dominance (SRHND) or sequen- 
tial dominance [8]. This additional constraint I2 = I3 will henceforth be applied in 
obtaining phenomenological solutions in the lepton sector. 

Apart from the obvious benefit of considering the neutrino sector, it will turn out 
that the neutrino sector will constrain the absolute values of the charges under the 
f/(l) family symmetry, (not the charge differences,) due to the Majorana nature of 
neutrinos. This is due to the relations between the charges imposed by the relevant 
GUT constraints, or the extended GUT constraints, eq. (fHTjl for the extended SU (5) 
solution of section 12.4.11 and eq. PTj) for the extended Pati-Salam solution of section 
12.4.21 For example the additional constraint I2 = I3 implies immediately 



I,. - I, 



0, 



(59) 



in the SU (5) type cases from EqlTTl 



Here we would like to study the cases for which large mixing angles in the at- 
mospheric sector and the neutrino sector can be explained naturally in terms of the 
parameters of the U{1) class of symmetries that we have constructed in the previous 
sections, under the framework of the type I see-saw mechanism together with the sce- 
nario of the single right handed neutrino dominance (SRHND). We refer the reader 
for a review of this scenario to [8]. Here we make a brief summary of the results 
and apply them to the present cases. In the type I see-saw the mass matrix of the 
low energy neutrinos is given by ^ v'^Y" M^'^Y^'^ , where Y" is the Dirac matrix 
for neutrinos and Mr is the Majorana matrix for right-handed neutrinos. If we have 
three right handed neutrinos. Mi, M2 and M3, then for the right handed neutrino 
mass, in terms of U{1) charges we have: 



M, 



RR 



^\ll+ni+hu\ 


^\li+n2+hu\ 


^\h+n3+hu 


^\l2+ni+hu\ 


^\l2+n2+hu\ 






^\h+n2+hu\ 




^\2ni+a\ 


^\ni+n2+o-\ 


^\ni+n3+a\ 


^\ni+n2+cr\ 


^\2n2+a\ 




^|ni+n3+(T| 


^|^^2+»^3+o■l 


£|2n3+<T| 



< s > 



(60) 



(61) 
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where the charges Ui are the U{1) charges of the right handed neutrinos, VRi and 
a is the U{1) charge of the field S giving Majorana masses to the right handed 
neutrinos. These charges are not constrained by the anomaly cancellation conditions 
Eq. (jni)-Eq. (fT^ of Section Q, at least in the SU{b) case, which gives some freedom 
in order to find appropriate solutions giving two large mixing angles and one small 
mixing angle for neutrinos. We expect S to be of order the scale at which the U{1) 
symmetry is broken, for example at Mp = Mpianck, or some other fundamental scale, 
such as the Grand Unification scale, Mq, for the solutions with an underlying GUT 
theory. 

Here we restrict ourselves to the cases in which Eq. 1)611) can be considered as 
diagonal, ~ diag{Mi, M2, M3}, for which we need in the (2, 3) block 

\n3 + n2 + (r\ > ■min{\2n3 + a\,\2n2 + a\}, 
2|n3 + ra2 + o-| > 12^3 + o-| + 12^2 + cr|. (62) 

The conditions in the (1, 2) block are analogous to the (2, 3) and also we need 

\ni + n-^ + a\ > max{\2n2 + cr|, |2?7,3 + a\}. (63) 

Now, there are two cases that we can consider here, which correspond to selecting 
which of the neutrinos will dominate. Mi or M3. For the later case the SRHND 
conditions are 



IM3I IM2I \Mi\ 

For the case in which Mi dominates we just have to interchange the indices 1 and 3 
in the neutrino Yukawa terms. 

For the case in which M3 dominates, at first order approximation, we have the 
following expressions for the neutrino mixings [8], 

^23 = ^, (65) 
-'33 

•yu (•Yv2 I V^2\ -yu (-yu •yu _ yu yu \ yu 

_ 121-^33 "'"-'23/ 13 1 33 32 -'22-'23/ ^ M2 /g^x 



V 

32 
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In terms of the Abelian charges the Yukawa elements are 

yu^^\U+n,+h^\^^\l'^+n,\^ = /. + Z,^ = _ 263, (68) 

where we have defined primed lepton doublet charges which absorb the Higgs charge, 
as shown. We can work here in terms of the primed charges, once they are fixed we can 
determine the original Abelian charges (unprimed). The approximation in Eq. ()fi7j) 
corresponds to the case in which we have enough suppression of the second term in 
the expression for tz/12. In Eq. (lUUj) the second term can be neglected sometimes, 
depending on the ratio M3/M2. The heaviest low energy neutrino masses are given 
by 

u2 2\V,+n,\ 2 u2 2\V,+n,\ 2 

™ _ 3 ™ _ 2 _^ (gg) 



M3 ' M2 ' 

where we have written e^K^+nsl ^ + and a^^^^li^+nal = (c23K,2 - S2zY^2f . 
Thus the ratio of the differences of the solar to atmospheric neutrino can be written 
as 



m,. 



M3ci3 (1^22 - >^32^^23) 



y2_ _ ^^M'-23V-'22 -Z2^'^-^i) _ _p2-P3 



M2 c?2 Yg + Yg 

where 



£p^-f^ (70) 



Pk = I2/2 + '^fcl ^ |2^fc + for k = 2, 3. (71) 



Note that is then defined such that 



m,, ^ — (72) 



6 SU{5) solutions satisfying the GST relation 

In this section we shall continue to focus on the case of SU{5), where the quark 

Yukawa matrices take the form of Eq ]501 where, motivated by large atmospheric 

neutrino mixing, we shall assume = from Eq. (jSIJjl The purpose of this section 

is to show how the GST relation can emerge from SU{5), by imposing additional 

constraints on the parameters. ^ 

^Note that results in Section and in Section [3 apply to both SU{5) type and extended SU{5) 
models, as discussed above. 
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6.1 The quark sector 



We have already seen that the GST relation can be achieved in the u sector, mainly 
by allowing the parameters s and s' to have different signs. In the down sector to 
satisfy GST we additionally require: 



\kd + r'. + s\ 



\kd + s'\ 



\kd + r'd + s'\ - \kd\ > \kd + r'^ + s\ + \kd 
W^ + kdl > \kd\. 



\kd + s\ 



(73) 



The first of these equations ensures the equality of the order of the elements (1,2) 
and (2, 1) of the matrix. The second equation ensures that the element (1, 1) is 
suppressed enough with respect to the contribution from the signalization of the (1,2) 
block. This last condition is usually satisfied whenever \kd + r'^ + s'\ > \kd + r'^ + s\ 
is satisfied. Finally the third condition ensures a small right-handed mixing for d- 
quarks and a small left-handed mixing for charged leptons. Now in order to satisfy 
the relations 



'12 



rrir 



'12 




A, 



(74) 



we need a structure of matrices, in terms of just one expansion parameter e = 0(A), 
such as 



for which we have 



... 1 



d 
'12 





" ... 




— 








_ ... 


_ 


23 ~ ^ 


' ■='13 ~ 





(75) 



(76) 



rric 4 rus 2 2 

vat fT^b 

in agreement with observed values for quark masses and mixings for e = \. 

Now we can proceed as in the example of Eq. where s' + s is fixed to be ±3. 
In this case we see that we can have plausible solutions in the up sector by allowing 
half integer solutions 



\s' + s\ = 13/2, 6, 11/2. 



(77) 
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We will refer to these solutions as Solution 1, 2 and 3 respectively. Note that only 
the charge differences are constrained here, the actual charges are not. 



Solution 1, \s + s'\ = 13/2, 



e 



35/2 pl3/2 ^35/4 



e 



e 



e 

13/2 ^9/2 
35/4 ^9/4 



r9/4 



^69/4 ^25/4 ^25/4 



^25/4 ^19/4 
pl7/2 p5/2 



,19/4 
r5/2 



for 



r', = h 

Solution 2, \s' + s| = 6, 



11, s 



1' ' 



35 
4"' 



— , kd 



-, or 



U = -11, s = s' 



35 i. 



-31/2 p-11/2 pll/2 



e 



e 



11/2 ^9/2 
15/2 ^5/2 



for 



10, 



Solution 3, |s + s' 

— 



^:. = ^i-^3 
"'^ = 11-1, 

= 11/2, 



^29/2 ^11/2 ^29/4 



-10, 



2, s' 



8, k, 



-8, fed 



,9/2 
,5/2 



5 

'2' 
5 

2' 



or 



^11/2 ^7/2 
^29/4 ^7/4 



^7/4 
1 



^31/2 ^21/4 ^21/4 
p21/4 ^15/4 ^15/4 



,33/4 



for 



/1-/3 
/1-/3 



41/4, s = s' 
7 



29 

T' 

29 



= -2, 



or 



4' 4 ' 



fcrf = 2. 



(78) 



(79) 



^0) 



^1) 



^2) 



(83) 



-41/4, s 

All the previous solutions Eqs. (f7^ -Eqs. (jHHj) lead to small tan/5 (0(1)), due to the 
choice of kd. To find solutions such that tan/5 is 0(10) is more difficult, due to the 
requirements in the up sector, but we have found the following solution 



/ / - 19 

'1 ~ '3 — TT' ^ 



r'd 



h-h 



-19 



-2. = 

2 

2, s' = -y, A;^ 



3 

■-, or 
2' 

3 



(84) 
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6.2 The neutrino sector 



Now we construct solutions for the lepton sector constrained by the requirements 
from the quark sector in the previous subsection, where we assumed = h — h = 0, 
and determined the charge differences r'^ = h — h that agree with the GST relation. 
Indeed it is convenient to label the solutions in the previous subsection by the value 
of r'^ = li — l2- Here we find the charges nj,/j, and a which satisfy the conditions 
arising from the neutrino sector, Eqs. (j65ll67|) . ^. In order to satisfy Eq. the 
most natural solution to achieve ^^'2 large is to have 

|/i +^2! = + ^2|. (85) 

The simplest solution is to assume that n2 = 0. Since l[ and are related through 
r'^ = li — I3 = l[ — Z2 the solutions to this equation are: 

r', = 0^ (86) 

I'l = | = -/2- (87) 

Since none of the solutions found in the previous subsection had r'^ = 0, we have to 
work with the second solution in Eq. (jS7j). However, we do not need to solve Eq. (jHTj) 
exactly, so we are going to perturb away from it, by keeping n2 7^ 0, but we expect it 
to be small in comparison with l[ = — /g- Then we write: 

P12 = \l[ + n2\ - \1'2 + n2\ (88) 

So is 0(£:^^^). The solution Eq. (jH7|) implies that l[ and I2 should have opposite 
sign, so we choose the case l[ > (the other case is similar). Since is large for 
all three GST solutions, and n2 should be small in order to satisfy Eq. (j88p . we can 
see that + ^^2! = — (^2 + ^^2), and \l[ + 712] = I'l + n2 for all the solutions from the 
previous subsection. Putting these relations into Eq. (|HHj) we get: 



n2 = — . (89) 

^The condition I2 = ^3 is a requirement of the class of see-saw models that we are looking for, 
single right-handed neutrino dominance (SRHND). Note that here we can also have I2 — — Z3 which 
then forces 713 = for ^ 0, in which case the solutions will be even more restricted. 
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So when we choose pu, n2 is determined. Now for the t^'s mixing, which should be at 
most 0(A), from Eq. (jMj) we need 

+^3] > + n3| ^ 723 > 0, (90) 

hence let us define P13 by: 

Pl3 = l^i+^sl - 1^2+^31, (91) 

We assume that the first term in Eq. ()66|) dominates. Then t^g ^ e^'^^/\/2. ^ By 
applying the same logic that led to Eq. (jH^ . we achieve: 



So fixing pi3 > 1 we fix 723. Now we need to impose the conditions under which 
we can have an appropriate value of Eq. ()70p. First note that in order to achieve 
m^3 = O(10-2)eV: 

for < S >= Mp, ^ 6 X 10"^ eV we need ~ 10^ 

< S > 

for < S >= Mg, ^ ^6x 10-^ eV we need ~ 10, (93) 

where has been defined in Eq. fj71|) . In terms of powers of A, we have A~^ — A^^ = 
O(IO^) - O(IO^) for < E >= Mp and A-\ A'^ = 0(10) for < S >= Mq. This 
corresponds to the following requirements: 

for < S >= Mp, p3 = (-4, -7) (94) 
for < S >= Mg, ^3 = (-1, -2). (95) 

We can conclude that for zero ^2, from Eq. (jU^ . since > 0, so must a be positive. 
Then we can write the power p2 — Ps {1711,2/171^.^ ~ ^pa-ps^ g^g follows: 

P2-P3 = -2(/; + 112) - {2112 + a) + (2^3 + a) T 2(/; + ^3). (96) 

The uncertainty in the final sign comes from whether IZ2I > If this is the case 
then we get: 

P3 -P2 = 4(n2 - n3). (97) 



^We have checked that this is indeed true for the solutions that we find for n2,n3 later in this 
section. 
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Otherwise we end up with 



P2-P3 = -4(/2 + ^2) (98) 

The second form is of no use to us, since we know that — is big for the models we 
are considering, and since n2 is small we can not get an acceptable mass ratio for iriy^ 
to m,^3. For the first form, Eq. (j^ . we need n2 ^ 0, because substituting Eq. 
into Eq. ()97|1 we have P2—P3 = 2pi3 — 4n2 and we need pis > 1 so for n2 = we have 
P2-P3> 2. 

With the above requirements then we can see that the parameters and n2 do 
not depend on r^. The only parameter which depends on this is a, through Eq. (|^. 
using the fact that I2 = —r'^/2. This also fixes the scale at which the U{1) should be 
broken. So, independently of r^, we have the following solutions 

1 ^ 3 11 

P12 = ^, Pi3 = 1, P2 - Ps = 2 ^ ^^2 = g , ^3 = 2 ; 

1 11,, 

P12 = 2 ' Pl3 = 1, P2 - P3 = 1 =^ n2 = 713 = -. (99) 

We can write the approximate expressions of mixings and masses in terms of the 
above results and the coefficients a^j of 0(1), 

a5'3£i2-3| _ a.^^|2n2| 

''12 



f — 

''23 ~ 


"23 
"33 


''13 




2 

^-23 


('^22 ~ 




r-v 2 

•-12 





V«33^ + «23^ ' (^23^22 " ^23032) ' 



^.2, ^^3 = 7^K3' + «23')^'^^'- (100) 

+ '^23 ) 

As we have seen above, the charges a are constrained by the differences r'^, the 
requirements of Eq. ()96|) and the solutions to Eq. (jU^ . which have the same value for 
n3, so for these two sets of solutions we have the same value for a. We write down 
these solutions for < S >= Mp in Table ^ and for < S >= Mq in Table (jD). 

The solutions presented here satisfy the conditions of the single neutrino right- 
handed dominance, Eq. (jU^ . which relate second and third families. For the first 
and second family we need similar conditions, which are safely satisfied whenever 
2ni > 2n2 > —cr for {2ni + a) positive. Thus rii is not completely determined but we 
can choose it to be a negative number between —o 12 and 0. 

Now that we have determined the conditions that the charges l\ and need to 
satisfy in order to produce SRHND solutions we can determine the Cj and U charges. 
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ool. 




^2 


^^3 


Pa 


a 


iW3 [GeVJ 


1 


11 


1 

8 


1 
2 


(-4,-7) 


(14,16) 


C(10 j, C(10 ) 


1 


1 1 

11 


1 
4 


1 
2 


(-4,-7) 


(14,16) 


C(10 ), C(10 j 


2 


10 


1 

8 


1 
2 


(-4,-7) 


(13,14) 


0(10^1), 0(10^) 


2 


10 


1 

4 


1 
2 


(-4,-7) 


(13,14) 


O(10ii),O(109) 


3 


41 
4 


1 

8 


1 
2 


(-4,-7) 


/53 57\ 
I 4 , 4 J 


O(10i0),O(108) 


3 


41 
4 


1 

4 


1 
2 


(-4,-7) 


/53 57^1 
I 4 , 4 J 


O(10i0),O(108) 


le 3: S at Mp 


for the solutions satisfying the GST relat 


Sol. 






^3 




O" 


M3 [GeV] 


1 


11 


1 
8 


1 

2 


(-1,-2) 


(10,11) 


0(10^) 


1 


11 


1 
4 


1 
2 


(-1,-2) 


(10,11) 


O(10«) 


2 


10 


1 
8 


1 
2 


(-1,-2) 


(10,11) 


O(109),O(10i0) 


2 


10 


1 
4 


1 
2 


(-1,-2) 


(10,11) 


O(109),O(10i0) 


3 


41 
4 


1 
8 


1 
2 


(-1,-2) 


/37 41 N 
V 4 ' 4 / 


O(109),O(108) 


3 


41 
4 


1 
4 


1 
2 
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Table 4: E at Mg for the solutions satisfying the GST relation. 



which are in agreement with the cancellation of anomalies, Eqs.(jS]|Hl), ^^*1 ^^^^ de- 
termines the matrices V^, and V^. In Section |3] we presented the conditions that 
the fermion mass matrices y", y^, and need to satisfy in order to give ac- 
ceptable predictions for mass ratios and mixings but without specifying the charges. 
The charges are then determined from r'^ and k^. We start by rewriting kd using the 
SU{b) charge relations, and the fact that = + hu- 



kd = q3 + d3 + hd = e3 + h- hu = 63 + (/g - h^) - K 



(101) 



Then we use the fact that A;„ = = 2e3 + h^i and we can solve for 63 in terms of kd 
and ( using Eq. (jHTj): 

2kd + 



es 



10 



(102) 



Once we have 63, and ^3, we can get ^3 since hu = — 2e3. From there, we can calculate 
the other charges from s,s',r,r' using Eq. pUj) and Eq. ()47|) . 

The charges calculated in this way are laid out in Table El 
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Table 5: Charged lepton charges for the SU{^) type solutions with m = u = 
satisfying the GST relation. The fits are discussed in section El 

6.3 Solutions for the extended SU{5) case with u + v ^ 

For this class of solutions, it is clear from Eq. and Eq. that the quark sector 
results will be unchanged. This happens since s,s',r,r' are blind to whether the 
family charges are related by the SU (5) relation, or the extended SU (5) relation, ku 
must always be zero, and the parameterization happens to leave k^ unchanged. Since 
ke is not unchanged, as discussed in section 1^.4. H we need to find ke in order to know 
the structure of the electron Yukawa matrix. 

It is helpful to rewrite ke and kd, from the form in Eq. (PT|) by using Eqs. (fT^ EZI) 
and k„ = 0: 



kd = /3 + 3e3 + M + -m, 
ke = /s + 3e3 + M + m. 



(103) 



where we have written m + f = m, as we will discuss in Section 110.31 m can be 
determined such a that the effects of the breaking of U{1) in the /i term are of order 
< m3/2. But on the other hand we need to keep the observed relation at low energies 
rrih = 0{mT-), so either m has to remain small or be negative to achieve \kd\ = 0{\ke\). 
In the present case the matrix has exactly the same form as in Eqs. ()39|) and 
has the form 



e 



e 



|s + fce| 



(104) 



With h = h, which determines the solutions of the charges Cj and k compatible 
with the condition Vd = h — h = 0, the discussion follows exactly as Section 
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because there we have not referred to other parameters than to kd, r, r', s and s' 
without specifying their relations with the charges cancelhng the anomahes. 

In this case, the analysis that leads to Eq. ()l()2j] must be repeated, but accounting 
for the fact that instead of the SU{5) relation between the charges, we must instead 
use the extended SU{5) relation between the charges. In this case, we find that: 

4 10 4 

kd = 3e3 + l3 + u + -{u + v)- 2K = Seg + ^3 + y « + 3^ (105) 

where we have used that l[ = li — 2e3 — u. / • is defined in such a way that U + rij + hu = 

r' 

l[ + rij. Using again the fact that l'^ = = find that: 

1 . , , 20 8 , , , 

e, = -{2k, + r',--u--v) (106) 

Using these results, and the values of s, s', rj, r^, we can find the charges in TablelHl 
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Table 6: Charged lepton charges for the extended SU (5) solutions with m = u + v = 
1/2 satisfying the GST relation. The fits are discussed in sectional 



7 SU{b) solutions not satisfying the GST relation 
7.1 The quark sector 

As we can see the GST relation puts a constraint on the opposite signs of s and s' 
and on the difference of = h — h. If we do not impose these requirements, allowing 
all the numbers s, s', r, r' and kd to have the same sign, positive or negative, we can 
factorize the k^ factor out of the matrix and so can write the down matrix in the 
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form 



r\kd\ 



^\s'+h-h\ ^\s'\ ^\s'\ 
^\s+li-l-i\ ^\s\ ^\s\ 



Ah-hl 



1 



1 



(107) 



In this case we do not have the restriction |s + Zi — /3I = |s'| so the parameter li — I2 
is not fixed by these conditions. In these cases is not constrained so it can acquire 
a value in the range ~ (0,3) for different values of tan/5. In these cases all positive 
or all negative charges, the cases which reproduce quark masses and mixings are for 



For Is I 



3 we have 
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(108) 



For Isl 



4 we have 
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(109) 



rlh-hl 1 1 



(110) 



From Eq. (jl09p and Eq. (jllO|) we can check if certain differences of leptonic charges can 
yield a suitable quark phenomenology. From Eq. ()5H1 we can see that in the cases of 
having all charges li and either positive or negative, then all the terms contributing 
to the first eigenvalue of V^, yi, will have the same power, as we mentioned earlier. 
So the difference here is constrained to reproduce an appropriate ratio nid/ms- Let 
us take here for definitiveness the case for positive charges (the negative charges case 
is completely analogous). Thus for s = 2, s' = 3, we have 



rrid 



(A2,A3/2) 



so in this case we have r'^ = 1, 3/2. For the case s = 2, s' = 4, we have 



(111) 



(A^A3/2) 



(112) 
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we do not want = as it will give somewhat large contribution from the (3, 1) 
element of the matrix to the eigenvalues. So for this case r'^ ~ 1/2. In this case 
we have the following matrices for Eq. 
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respectively for 



3/2. For Eq. (fTTT]|l we have 

.9/2 



e 

£5/2 

el/2 



e2 e2 



(113) 



(114) 



for = 1/2. These solutions work for kd G (0,3), depending on the value of tan/5, 
these matrices yield acceptable phenomenology in both charged the lepton sector and 
d quark sector. 



7.2 The neutrino sector 

As we have seen in Section [7[TJ in these cases is constrained to be G (1,3/2) 
for (s, s') = (2,3) and ~ 1/2 for (s, s') = (2,4) but let us leave it unspecified for 
the moment. We consider here the case of all the parameters related to U and Cj 
positive. In this case we require that all the neutrino charges, rij to be negative but 
a positive. We proceed as in Section ()fi.l|l . in order to identify the charges rii and 
a. In principle we need = 5K2+"2| i^^^^ require /^,/2 > so now the 

appropriate solution to this would be 

/'i = r:„/^ = 0, ri2 = ^. (115) 

However in this case, as in the case of section ()(i.l|l . we will only be able to produce 
m^j/myg ~ So we work with a solution of the form Eq. (jHHj) . For this case we 
then have 

l[ = r'd, I', = 0, = ^i^. (116) 

Note that in this case the charges U are positive because h = kd — 3e3 and here 
63 = kd- For ti3 we also make use of the parameterization of Eq. ()9H). Assuming that 
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\r'J > {risl 



ns = (117) 



In order to achieve an appropriate ratio for 171^^/171^.^ we need now the conditions 
2^3 + cr > 0, 2^2 + cr > 0, Z2 + < 0, Z2 + ^3 < 0, for one of the last two inequahties 
the equahty can be satisfied, but not for both. For this case, we have also P2 — Ps = 
4(^3 — n2) and using Eq. pi6|) and Eq. ()117|) we have P2—P3 = 2(pi3 —pu)- We can 
also choose the parameters pu, P13 and p2 — Ps as in Eq. but now and n2 are 
given by Eq. pi6p and Eq. piTj) . Thus we have 

P12 = -, Pl3 = 1, P2-P3 = ^ 

^^2 = ^- ^<0, n3 = ^- ^<0^r:,>l, (118) 

P12 = ^, Pl3 = 1, P2 -P3 = 1 

^n2 = i- ^<0, n3 = ^- ^<0^r^>l. (119) 

In Section (jH) we determined the approximate values for r^. For (s, s') = (2, 3) we 
can have = 1,3/2 while for (s, s') = (2,4) we have r'^ ~ 1/2, which however is 
not in agreement with the conditions of Eq. pi8j) and Eq. ()119j) . The approximate 
expressions of mixings and masses in terms of the above results and the coefficients 
a^j of 0(1) are as in Eq. (llOUp . except for and t'(2 which now read 

^ ekd+"3|-|rt3| a^,£:(l'"d+"2l-l"2|) 



^13 - / „ 9 . 770 ' ''12 



^/^^7T^ ' (C23a^2 - ^230^2) 



(120) 



We have listed the possible solutions of Table (0) for Eq. piSj) at (S) = Mp and in 
Table © for (S) = Mg. 



7.3 Solutions for the extended SU{5) case with u + v^O 

We do not present any charges for this class of solutions, but here is how the charges 
would be calculated. In this case, the analysis is carried out in the same way as in 
section inini The only subtlety is that the relation linking I2 to is different. Instead, 
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Table 7: S at Mp for the solutions not satisfying the GST relation. 
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Table 8: S at Mq for the solutions not satisfying the GST relation. 

we have, from Eq. ()116p that ^2 = ^3 = 0. Putting this result into Eq. ()1U5|) . we 
achieve: 

1 . , 4 8 , , , 

63 = Y5(2/c, + -w--t;). (121) 

From 63 and Z3, the other charges may be calculated using the known values of 
s, s', Trf, r[^, M andf , by using the extended SU{^) charge relations, Eq. (jHTj) and the 
simplified parametrization, Eq. PHj) . 



8 The non-SUip) Cases 



8.1 Solutions for w = = 0, in the Pati-Salam case 



With I2 = h, in this case we have s 
matrices are 



^2 — = h — h = then the charges of the 



dA\ 



h-h + di- 63 /l - /3 + 62 - 63 h - I3 

ei - 63 62-63 

61 - 63 62-63 

li + l3 + ei + 63 h + l3 + e2 + 63 /i + /3 + 2e3 

2/3 + 61 + 63 2/3 + 62 + 63 2/3 + 263 

2/3 + 61 + 63 2/3 + 62 + 63 2/3 + 263 



(122) 
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Table 9: Charged lepton U{l)x charges for the solutions m = f = not satisfying the 
GST relation. The fits are discussed in section IHl 

In this case the U{l)x symmetry does not give an appropriate description of fermion 
masses and mixings, however it can be combined with non-renormalizable operators 
of the Pati-Salam group, [30], in order to give a good description of the fermion 
phenomenology. 



8.2 Solutions for u + v = {) 



One trivial example of non- SU (5) cases was given in section 1)2.3.11) for the solution 
u + V = 0. We proceed as in the section ( 0])- in order to analyze the appropriate 
phenomenology. We are interested in the cases I2 = h, this together with the condition 
of 0(1) top Yukawa coupling give us the following matrices of charges, which are 
derived with the appropriate substitutions in Eq. (jHUl-Eq. 



h + ei 

^ + 62 - 63 



+ 61-63 



/3-/l 



h + 61 

+ 63- 

h-h 



63 



62 - 63 



3 h 63 — Ci 

63 - 62 





62 - 63 
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3 +^3 
63 - 62 





61 



h + 61 h + 62 h + 63 

62 - 63 62 - 63 
62 - 63 62 - 63 



(123) 



Due to the form of the charges in the up and down quark matrices, first at all we 
would need two expansion parameters: and e*^. But with this structure alone it is 
not possible to account simultaneously for appropriate mass ratios of the second to 
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third family of quarks and for an appropriate Vcb mixing. So in this case just with a 
U{1) it is not possible to explain fermion masses and mixings in the context of the 
single neutrino right-handed dominance, SNRHD. 

9 Numerical fits of masses and mixings 
9.1 Fitted examples 

In this section we present numerical fits to some of the examples detailed in Sections 
()(jp7|l and we compare the results with a fit of a generic SU{3)-like case [6]. The 
simplest way to construct the lepton Yukawa matrices from the charges is to first 
calculate h^^d- We extract from k^., I3 and 63 from ke = I3 + + h^. In general, 
we can use Eq. fll03|) to obtain: 

hu + hd = m = 3{kd - ke) (124) 

This is then enough to construct the lepton Yukawas from the appropriate line of 
the table (jSl or IH)) of the lepton and Yukawa family charges. Below we specify the 
examples that we have chosen to fit. 

Fit 1: SU{5) type solution (u = v = 0): example satisfying the GST relation 

This takes GST solution 2, (Eq. (IHUj) ) in the SU{5) type case, with u = v = 0. The 
charges Cj, and n2,3 are read off from the fourth line of Table In principle, the 
value of a would be read off from either Tabled] ( for neutrino masses generated at the 
GUT scale ) or Table El ( for neutrino masses geneated at the Planck scale). However, 
these tables allow for a range of a; for this fit, we take a = 21/2 for GUT scale 
neturino mass generation, and cr = 29/2 for Planck scale neutrino mass generation. 

Then, up to a and rii : —a/2 < ni < 0, the Yukawa and Majorana matrices are: 
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(125) 



Fit 2: Extended SU{5) solution {u + v j^O) satisfying the GST relation 

This takes GST solution 2, (Eq. (jHUjl ). in the extended SU{5) case with u + v 0. 
The charges /j, Cj and 71,2,3 are read off from the third line of Table El The values 
of cr taken are a = 19/2, a = 29/2 for GUT scale and Planck scale neutrino mass 
generation respectively. Again, rii is taken to lie in the region —ajl < rii < 0. ®. 
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Fit 3: Extended SU{5) solution {u + v j^O), satisfying the GST relation 

This takes GST solution 3, ( Eq. (jEHI)), in the extended SU{5) case with m + t> 7^ 0. 
The charges /j, Cj and n2,3 are read off from the fifth line of table IHl The values of a 

^The difference between Fit 1 and Fit 2 is that the charges (Tables (0) and ^ respectively) are 
determined in a different way and hence the value of the effective parameter expansion e is different 
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taken are a = 39/4, a = 55/4 for GUT and Planck scale neutrino mass generation 
respectively, ni lies in the region —a/2 < rii < 0. 
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Fit 4: SU{h) {u 



0) solution not satisfying the GST relation 



Here we present a solution non satisfying the GST relation of the form of Eq. p()9j) 
for /i — Z3 = 1, which corresponds to the set of charges of the first line of Table 0- 
We also fix here the expansion parameter e = 0.19, using the FI term. The high 
energy Yukawa and Majorana matrices are: 
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Fit 5: S't/(5) (m = f = 0) solution not satisfying the GST relation 

Here we present another solution non satisfying the GST relation of the form of 
Eq. ()lU9p for li — 13 = 3/2, which corresponds to the set of charges of the second line 
of Table We also fix here the expansion parameter e = 0.185, using the FI term. 
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The high energy Yukawa and Majorana matrices are: 
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9.2 Details of the fitting method 

One of the purposes of these fits is to compare which solution fits the data best 
while constraining the abritary coefficients to remain at 0(1). We therefore choose a 
minimization routine to find these 0(1) coefficients and compare the numerical values 
for the different solutions. In the quark sector we use eight experimental inputs in 
order to determine the parameters (coefficients or phases): 

Vub/Vcb, Vtd/Vts, Vus, Im{J}, mu/mc, rric/mt, md/rris, rris/mb. (130) 

We explain in the Appendix how this fit is performed, the important point is that 
we can only fit eight parameters and the rest need to be fixed. The minimization 
algorithm has been optimized to fit the solutions satisfying the GST relation because 
the number of parameters is close to eight. We also fit examples of the non GST 
solutions but since there are more free parameters in this cases (mainly phases) it is 
un-practical to make a fit by fixing so many free parameters. So we present particular 
examples in these cases which do not necessarily correspond to the best x^- 

In the lepton sector we perform two fits, one for the coefficients of the charged 
lepton mass matrix and the other for the coefficients of the neutrino mass matrix. 
We do not perform a combined fit for the coefficients of and because the un- 
certainties in these sectors are quite different. While the uncertainties in the masses 
of the charged leptons is very small, the uncertainties in lepton mixings and quan- 
tities related to neutrino masses are still large, such that we cannot determine the 
parameters involved to a very good accuracy. 
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The quantities used for the fit of the coefficients of the charged lepton mass matrix 

are 

^, ^, (131) 
rur 

such that we can just determine two parameters, a\2 and 022, but for the cases pre- 
sented here this is enough. In order to do the fit for the coefficients of the neutrino 
mass matrix we use the observables 

/ / / I ^^^Sol I / N 

^235 ^13' '^125 I '\i '^1^3 (13^) 

I '^atm I 

where we relate ^33 to the atmospheric mixing, to the solar mixing and t'^g to the 
reactor mixing. In this case we are going to be able to fit just five parameters. For 
this reason and because the uncertainties in the above observables are significantly 
bigger than the uncertainties in the quark sector, the fits of the coefficients of the 
neutrino mass matrix have large errors and they may leave a room for other solutions 
once the experimental uncertainties improve. Since we only have an upper bound for 
the reactor angle, t\^, we fit the solutions in the neighborhood of this upper bound. 



9.3 Results of the fits 

9.3.1 Fit 1: SU{b) {u = v = 0) example satisfying the GST relation 

This is a SU{^) type solution, and hence m = f = 0, which satisfies the GST relation. 
The textures are as laid out in Eq. (jl25p . 



Quark sector 

We can use the expressions Eq. flFIj) and Eq. adapted to the solution of Eq. (jST|) 
in order to fit the Yukawa coefficients, along with the appropriate phases entering 
into the expressions of mixings. The expansion parameter e is determined with the 
Fayet-Iliopoulos term and the appropriate charges cancelling the anomalies, for this 
case its value is e = 0.183. The parameters that we fit are the real parameters 

^121 ^21,1 ^221 ^121 ^11,1 '^23' '^32; C0S($2), (133) 

which enter in the expressions of mixings and masses, Eq. (j^T|) -Eq. ()54|) . Note that 
in these expressions the coefficients af • can be complex but for the fit we choose them 
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real and write down explicitly the phases. We are free to choose the parameters 
to fit. However we need to check which are the most relevant parameters to test 
the symmetry. Thus we follow this as a guideline to choose the parameters to fit 
and leave other parameters fixed. Due to the form of Eq. (jHT|) the mixing angles in 
the (2, 3) sector of both matrices contribute at the same order in the Vckm matrix 
mixing, s^g = [023 — 0236**^23 [e^, so we have decided to put a phase here. In the ^"2 
diagonalization angle and the second eigenvalue of the combination 022^**^ — (223^ 
appears, so we have chosen as well to include a phase difference there. The fixed 
parameters are then 



^22' 



$1, $3 



^23 ' 



(134) 



where $1 has the form of Eq. ()56|) and the phases $3 and $^23 can be written as ^ 

$3 = 022-20L <^'X.3=(033-023)-(033-</>23)- (135) 



The results of the fit in the quark sector appear in the second column of Table 
fllO|) . Given these results we can think that the structure of Yukawa matrices has the 
following form 



* 


l/i2e'*i 


1/13 


Z/i2e'*^ 


l/22e'*^ 


1/23 


1/13 


1/23 


1 



yd 



l/2ie**2 



1/126 
1/22 
1/32 



j*2 



Vise 

1/23 

1 



(136) 



where yij denote real elements and we have associated the phases to particular 
elements of the matrices. Note that we need three phases to determine the amount 
of CP violation experimentally required because in all the fits we found $^23 = 0- If 
this phase was not zero then it could have been associated to the Y23 element. The 
entries marked with * cannot be determined because they are not restricted by masses 
and mixings, due to the structure of the Yukawa matrices. The value of 1/21 e 
determined indirectly because we need to satisfy the GST relation so 
both up and quark sectors. 



IS 



ti2 for 



Lepton sector 

We have fixed the coefficients of in the quark sector and now we can use the results 
for the charged lepton matrix F^. The masses of the charged lepton are obtained 
^In terms of the Pi phases appearing in the diagonaUzation matrices, Eq. 1176(1 . we have $1 = 
^, $2 = -f3l ^ and <^x,, = {[3^ ^ - ^) - {P^ ^ - /JJ' ^). 
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Quark Fitted Parameters 


GST sol. 2 


GST sol. 2, M,f 7^ 


GS i sol. 3, M, f 7^ 


Parameter 


BFP Value 


BFP Value 


BFP Value 




2.74 ± 0.61 


1.04 ± 0.19 


2.74 ± 0.71 


^23 


1.68 ± 0.17 


1.34 ± 0.13 


1.41 ± 0.18 


^22 


1.08 ±0.18 


1.05 ± 0.11 


0.70 ± 0.23 


<2 


0.93 ±0.15 


0.55 ±0.20 


0.74 ±0.13 


<3 


0.29 ±0.21 


0.30 ±0.14 


0.74 ±0.17 


'^23 


0.79 ±0.10 


0.70 ±0.13 


0.66 ±0.35 


"32 


0.48 ±0.17 


1.28 ±0.32 


1.28 ±0.58 


C0S($2) 


0.454 ±0.041 


0.456 ± 0.041 


0.547 ±0.424 


Quark Fixed Parameters 


e 


0.183 


0.217 


0.154 


'^22 


1 


1 


1.4 


C0S($3) 


0.8 


0.83 


0.8 


C0S($X23) 


1 


1 


1 


$1 


7r/2 


7r/2 


7r/2 








1.47 


2.41 


4.32 



Table 10: Quark fitted parameters for the examples of Section|HJ. The second colmim corresponds 
to the Solution 2 in the SU{5) {u = v ^ 0) case, the third column to the Solution 2 in the u ^ —v ^ 
case. The fourth column presents the fit to the Solution 3 in the u ^ —v ^ case. 

through the SU{5) relations, ensuring the correct value of charged lepton masses, 
once the masses of the d-quarks are in agreement with experimental information. 
Thus in this case we perform a fit just for coefficients of the neutrino mass matrix, 
y^, using the ratio of neutrino mass differences (solar to atmospheric), the mass of 
the heaviest neutrino and the lepton mixings, which have a contributions from both 
the charged leptons and the neutrinos. Here the relevant parameter that we need 
from the quark sector is because the tangent of the angle diagonalizing on the 
left is related to this parameter: t^^ = (123 (^32- Since this is an 0(1) mixing we 
have to take it into account for the results of the Umns mixings, thus we have 

|„e „u „—i(pXn^ a'^ n'^ \ 

_ l'^23'^23^ ■^23'"23l (137) 

23 q<^ _|_ r'^ /-e piipXoi I ' 

l''23''23 ^ '-23'-23'^ I 
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where we use the expression Eq. to determine S23 ci-nd C23, and the approximation 
^23 — ^321 'Pxn is a phase relating e and z/ mixings in the (2, 3) sector [38]. We denote 
the Umns angles by the superscript / and by e and v the charged lepton and neutrino 
mixings respectively. The mixings ^'^^3 and ^^2 essentially given by the neutrino 
mixings, Eqs. pOO|) . so we fit these mixings according to Eq. (jUBj) and Eq. (p7j) 
respectively. We note from Table (fTTj) that in the lepton sector we need two phases, 
0X23 ^'^d (j)^ . The phase 0^23 can be associated to the charged lepton sector and we 
can put it in the 1^23 entry. The second phase, cj)^ can be assigned to ^22- We fit 
the mass ratio and the heaviest neutrino state using their expressions appearing in 
Eqs. fjlOOp . The results for this fit appear in the second column of Table (fTTj) . 



Neutrino Fitted Parameters 


Parameter GST sol. 2 


GST sol. 


2, M, w 7^ 


GST sol. 3, u,v ^ 




Mp 


Mg 


Mp 


Mg 


Mp 


Mg 




BFP value 


BFP value 


BFP value 


BFP value 


BFP value 


BFP value 


'^23 


0.75 ±0.79 


0.67 ±0.61 


0.21 ±0.25 


0.85 ±0.27 


0.30 ±0.18 


0.40 ±0.15 


«r3 


1.41 ± 1.32 


1.36 ± 1.10 


0.97 ±0.47 


1.25 ±0.63 


1.02 ±0.50 


1.45 ±0.70 


"12 


2.23 ±0.92 


2.10 ±0.81 


1.25 ±0.29 


2.08 ±0.69 


1.35 ±0.34 


1.97 ±0.45 


'^22 


1.84 ± 1.37 


1.96 ± 1.92 


1.23 ± 1.41 


1.98 ±0.79 


1.48 ± 1.31 


2.26 ± 1.6 


"32 


1.47 ± 1.93 


0.98 ±0.91 


0.65 ±0.70 


1.53 ±0.75 


0.53 ±0.78 


0.56 ±0.98 


Neutrino Fixed Parameters 


e 


0.183 




0.217 




0.154 




^22, 


0^2 = 0.48 




-1.6 




1.2 




'^33 


1 




1 




0.7 


1 


a 


29/2 


21/2 


29/2 


19/2 


55/4 


39/4 


C(0X23) 


0.29 




0.29 




1 


0.5 


cm 


-1 




-0.5 




0.86 


1 


n2,n3) 


(1/4, 1/2) 


(1/8, 1/2) 






0.44 


0.12 


1.67 


0.49 


2.16 


0.72 



Table 11: Neutrino fitted parameters for the examples of Sectional The second column corresponds 
to the Solution 2 in the SU{5) {u ^ v ~ 0) case, the third column to the Solution 2 in the u + v ^ 
case. The fourth column presents the fit to the Solution 3 in the u + v ^ Q case. Here c(y) is the 
cosine of the respective parameter. 
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9.3.2 Fit 2 and Fit 3: Extended SU{5) solutions with u + v satisfying 
the GST relation 

These are both extended SU{5) solutions, with u + v ^ 0, satisfying the GST relation. 
Fit 2 corresponds to the textures laid out in Eq. p26|) . and Fit 3 corresponds to the 
textures laid out in Eq. p27|) . 

Quark sector 

This section is completely analogous to the previous one, the only difference is in 
the value of e. We present here two examples. The first example corresponds to the 
first solution of Eq. (jHTj) . which we called Solution 2, and corresponds to e = 0.217 
according to the charges of the third row of Eq. ©. The second example corresponds 
to the first solution of Eq. (jH^ . which has been called Solution 3 and corresponds 
to £ = 0.154, according to the charges of the fourth row of Eq. 0. The fitted and 
fixed parameters are also those of the previous example, Eq. p33|l and Eq. p34j) 
respectively. The results for the quark fitting are presented in the third and fourth 
column of Table (fTUj) . respectively, so we can compare directly with the previous case. 

Lepton sector 

This case is different from the Section 19.11 because now we do not have the SU (5) 
relations. Instead the parameter is different from fc^, as explained in Section (j3)), 
and hence ^ {y^y . In this case we perform two fits, one for the coefficients of 
the charged lepton mass matrix, and another for the coefficients of the neutrino 
mass matrix, Y^ . 

For the Solution 2, taking into account the value of the charges, the second row 
of Table ©, and that m = u + v = 1/2 we have k^. = —8/3. We note in this case 
that since we need ^ m^, which are given by 

rrib = mtel*^'*!, ka = h + Ses + u + A{u + v)/3 

TTLr = mte^''^^, ke = h + 3e3 + u+ (u + v), (138) 

we expect the sum {u + v) to remain small. 

Now the coefficients a?,n and aio are not related but we can fix af,n in the neutrino 
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sector such that it is in agreement with the results from neutrino oscillation. We have 
performed a fit using the experimental information of the parameters of Eq. ()132|) . 
Here we have also used the expression Eq. ()137|1 in order to fit the atmospheric 
angle, the expressions Eq. and Eq. (jHTj) to fit t[.^ and (reactor and solar angle 
respectively) and the mass ratio and the heaviest neutrino state using their expressions 
appearing in Eqs. ()100|) . The results for this fit appear in the third column of Table 

m- 

Once the parameter 0^3 has been fixed we fit the parameters of the charged lepton 
mass matrix, of the form Eq. ()l()4j) and the other parameters as in the first solution 
of Eq. ()8H). In this case the relevant parameters and 022- However if we just 

fit the expressions 

I e 1 2 

_ 4/3 _ „e 2 



I / e e e \ '"^ 

^/i K'^22 ~ '^23'^32J 



_ _ '12 5 

I K"22 "23"32y'l 



«22 - "2342)^ . (139) 



the coefficients 0^2 and not quite 0(1) so we have to make use of a coefficient, 

c such that (032 — ct23'^32) ('^22 '^23'^32)/C' ^'S- c = 3, in order to have acceptable 
values for charged lepton masses. This fit is presented in the second column of Table 
fjl2|) . In this case the extra-coefficient needed for the fit is not really justified in the 
context of just a single U{1) symmetry. 

For the Solution 3, we have m = 1/2, = —13/6, according to the charges of 
the third row of Table (jHI). The fit of the coefficients of the neutrino mass matrix are 
completely analogous for Solution 2 and they appear in the third column of Table 
flll|) . The relevant parameters for the charged lepton sector are 

I "12 I 29/6 _ „e 2 



^ ~ "^12 ) 



me 

I ('^22 ~ '^23*^32) 

«22 - a23«32)^'/'- (140) 



For this case there is no need to invoke another coefficient as for the Solution 2. 0(1) 
coefficients in this case can account for the masses and mixings in the leptonic sector. 
Once the coefficient 0^3 is fitted in the charged lepton sector then we need to use this 
parameter as a fixed parameter in the fit for the neutrino sector but in this case the 
fit is not as good as for the previous solution. The results are presented in the third 
column of Table (fT^ . 
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Charged lepton Fitted Parameters 


GST sol. 


2, ^ 


GST sol. 3, M, 7^ 


Parameter 


BFP Value 


BFP Value 




0.56 ±0.006 


2.88 ±0.032 


'^22 


0.92 ±0.013 


1.87 ±0.013 


Charged lepton Fixed Parameters 


e 


0.217 


0.154 


'^23 


-1.6 


1.2 


'^32 


1.8 


1.2 






0.05 


1.2 X 10-5 



Table 12: Charged lepton fitted parameters for the examples of Section |S1 The second column 
corresponds to the Solution 2 in the u + v ^ Q case. The fourth column presents the fit to the 
Solution 3 in the u + v ^ Q case. 

9.3.3 Fit 4: 5'f/(5) type {u = v = Q) solution not satisfying the GST 
relation 

This is a SU{b) type solution, hence u = f = 0, which doesn't satisfy the GST 
relation. The charges are as laid out in Eq. ()128j) . 

Quark sector 

Here we also use the expressions Eq. fl^ and Eq. adapted to the solution 
Eq. ()109|) for r'^ = li — = 1 and check the fit with an exact numerical solution, 
which agrees with the fit to Eq. (|3T|) and Eq. within a 5% error. Since the fit 
can just fit eight parameters, in this case it is not possible to select out "the best fit", 
according to the criteria that we have used for the previous fits, so we present the 
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following solution for the coefficients of the up and down Yukawa matrices: 

0.58e-^'^/2 

0.43e-^'^/2 
e 0.8 



0.42 

0.58e-^''/2 

0.51 

-jO.5 



-iO.79 



0.86e-'i-2 

Q^4g-i0.5 



0.51 

0.43e-^^/2 

1 

Q 29e^o.48 ■ 
0.55e-*-^ 

„-i3.05 



(141) 



For this fit we have = 2.31. 
Lepton sector 

In the lepton sector, once we have done the fit to the quark masses, the SU{5) relations 
produce acceptable values for the charged lepton masses, what we need to care about 
are the mixings for the neutrino sector. According to the expressions for the mixings 

a' 



in the (1,2) and (1,3) neutrino sector, Eq. p20|) . now f^. 

tl2 = ai2^^^^/(c23a22 - 43«32))' (^2, ?^3 

mixings in the charged lepton sector go as t 



13 



"13' 



,u 2 

'33 



+ 



e 

12 



(—3/8,0). On the other hand the 
l«2i + 3a23«3i/«33k/3|a22 + 3a: 



and t^3 



ai^e/\a 



33 + I '^32 1 



32*^231 
I 



SO here these contributions are important to the Umns Su 



and mixings, identified respectively to the solar and reactor mixings, for example 
for we have 



'13 



I''12''13''13 



C^3(e^(^i-^i)r- 



'^23l'^12'^23"^13 + ^ 



e e \ 
12*23 J 



if}: 



^■^12^13 ^12813)323) 



(142) 



The mixing S23 is driven by the neutrino mixing s^. 



J J 

'23'-13 



''23''12''23 



23 



■^23^13*^23 I 



(143) 



Despite all the contributions to the mixings and s[2 we can reproduce the observed 
masses and mixings in the neutrino sector with 0(1) coefficients and with out any 
phase in this sector, we just use the phases of the right handed quark matrix, which 
are given by 



/5| 



ArcTan 



^32 
^22 



sin(( 




_cos((/)^3) 


+ \42? 


d \ j_ odR 
33J "T Pi 




d \ odR 
21) P2 





(144) 



45 



and are specified in Eq. ()14H) . The results of tliis fit are given in the second row of 
in Table 



Neutrino Fitted Parameters 


Non GST sol. 1 


Non GST sol. 2 


Parameter 


BFP Value 


BFP Value 


^23 


1.6 ±0.8 


2 ±0.9 


ah 


1.4 ±0.7 


0.9 ±0.3 


0-12 


1 ±0.6 


1.6 ±0.3 


^22 


0.67 ±0.27 


0.5 ±0.4 


Neutrino Fixed Parameters 


e 


0.19 


0.185 


0-23 
^33 
'^22 

a 


-3af2 = -1.2 

1 

2 

(4.5,0.5) 
(-3/8,0) 


-3a;f2 = -1-25 

1 

2 

(5,1) 

(-5/8,-1/4) 






(5.09,4.77) 


(4.78,3.79) 



Table 13: Neutrino fitted parameters for two of the non GST examples of Section [T] The second 
and third columns correspond respectively to solution 1 and 2 in the non GST SU (5) {u = v — 0) 
cases, for the first one we have used = 1 and for the second r'^ = 3/2. While we have fitted 
in the first case to saturate its current upper limit, we have allowed for the second case to be 
smaller than it. The first entry for a corresponds to the fit using Mp and the second entry using 
Ma] analogously for ■ 

9.3.4 Fit 5: SU{5) type (u = v = 0) solution not satisfying the GST 
relation 

This is a SU (5) type solution, and hence u = v = which doesn't satisfy the GST 
relation. The textures are as laid out in Eq. ()129p . 
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Quark sector 



Here we present the following solution for the case Vd = h — h 
coefficients of the up and down Yukawa matrices: 



|, in this case the 



0.5 0.6e-*^/2 

0.5 0.43e-*'^/2 
1 

1.82e-^2.28 



0.72 

0.76e-*i-i2 



-jl.57 



0.4e 



-i0.41 



0.5 

0.43e-^'^/2 

1 

0_29e*o-49 
0.55e-*°-^^ 

„-i2.951 



(145) 



For this fit we have = 2.10. 



Lepton sector 

The analysis of this fit is completely analogous to the Fit 4, the results of the fitting 
procedure is presented in the second column of Table [121 

9.3.5 Top and bottom masses and tan/3 

For these cases tan (3 and prediction, once the coefficient 033 is fixed through 

the value of rrit, rrit = Y^^v/ \/2. The values of 033, 033 and tan (3 for the cases presented 
in this section are given in Table (fT^ . We can see that for a natural value of 033 = 1 
we have acceptable values for tan (3 (which should be > 2) and af^ in any of the cases 
presented. 



9.4 Comparison to the 5f/(3) case 

In this section we present the comparison to a generic SU{3) case [6]. What we fit 
are the 0(1) coefficients of a Yukawa matrices of the form 







F^ 


F^ ' 


Yf = 




4 






F^ 


4 


1 
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tan p, 


033 and 033 




Parameter 


LrO 1 sol. Z, M= V = U 


Lto i sol. z, 7= U 


KjD I sol. C), f 7= U 


•^33 


^^ 1 Q/i ^ 


(1.1.3) 


(1,1.3) 


nd 
"33 




l'-'-^^+1.84? -'-•"^+0.08J 


(0 90-0.68 -r «9-0.09\ 
I'-'-^'-'+l.TO' -'-•"^+0.10/' 


tan p 




(3.00;?:^?, 1.07;°:°?) 


(3.oo;?f , 1.07;°:°?) 






(0.217, 5/2) 


(0.154, 2) 


Parameter 


iNonvjroi sol. i, — f — u 


Non GST sol. 2,u=v = Q 


^33 


(1,1.2) 


(1,1.2) 




•^33 


(c) 1 9-0.44 1 1 -0.03\ 

|^Z.iZ_,_Q gg, i.i^Q QgJ 


/9 -| 1 -0.45 -| q-0.06\ 
l^Z.±±_,_Q i.O_,_Q qq; 




tan P 


/q-o.ee 1 o-o.i\ 

l'J+1.61' -'-•"5+0.1/' 


/q-0.78 -1 9-0.2\ 
W+1.32) ^■^+0.2) 






) 


(0.19, 2) 


(0.185, 2) 





Table 14: Value of afg and tan/3 for the different models presented, once 033 is fixed using m*. 

where we allow two different expansion paramaters e„ and Ed and complex phases to 
reproduce the CP violation phase. It is enough to consider one different phase in each 
of the y" and matrices. Here we put the phases on Yf^ and y,^^ [20] , but we have the 
freedom to use other choices. We have used here as well the method of minimization 
that we have used for the C/(l) cases. The results of these fits are consistent with 
previous determination of these parameters, [19,20], taking into account the change 
induced by the change of the value used here for the parameter rric/ms — 15.5±3.7 and 
the different methods used for the determination of coefficients^^ . The fits presented 
here are the fits with the lowest possible because of the minimization procedure. 
We have not included here for the SU{S) case a fit in the neutrino sector because 
in the 5'C/(3)-like cases the neutrino sector requires more assumptions than in the 
analogous U{1) cases. 

Another important difference between the SU{3) and the U{1) cases presented 
here is that in the first one there are two parameter expansions and s*^ which 
have been fitted while in the U{1) cases there is only one expansion parameter which 
can be fixed by relating the U{1) symmetry to the cancellation of anomalies and the 
Fayet-Iliopoulos term. This has allowed that more 0(1) coefficients have been able 
to be fitted. 

lojn [19,20] rric/ms = 9.5 ± 1.7. 
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Quark fitted Parameters, SU(3)-like case 


Parameter 


BFP Value ±a 


BFP Value ±cr 


a'u22 


1.11 ± 0.55 


1.11 ± 0.07 




0.66 ± 0.32 


2.45 ± 0.20 




0.10 ± 0.12 


0.91 ± 0.15 


'^22 


0.74 ±0.10 


1.77 ±0.09 


'^23 


0.45 ±0.29 


1.18 ±0.12 


e" 


0.05 ±0.007 


0.05 ±0.007 




0.25 ±0.03 


0.16 ±0.02 


C0S(<I>2) 


0.516 ±0.1 


0.450 ± 0.045 


Quark Fixed Parameters, SU(3)-like case 




-1.25 ^= -0.87r/2 


1.120 ^ 0.77r/2 






0.972 


0.974 



Table 15: Fitted and fixed parameters for the 5f7(3)-like case. 

By comparing Tables ()10|1 and ()15|1 we can see that according to the minimization 
procedure and the criteria of 0(1) coefficients, the second case of the SU{3) solution 
fits better the data. However the U{1) solutions also have a good fit and taking 
into account the fact that for the neutrino sector we just have added the SRHND 
conditions, the fits in both of the U{1) cases presented are good. We can therefore 
consider that U{1) symmetries are still an appealing description of the fermion masses 
and mixings observed. Note that although the Solution 3 in the u ^ v does not 
fit the data as well as the Solution 2 (in either case, u = v = or not) in the quark 
sector, it does reproduce masses and mixings in the charged lepton sector. We have 
for this case ^ {Y'^)'^ but we have ^ m,- without introducing ad- hoc 0(1) 
coefficients in order to reproduce the appropriate mixings. 

Given the results of these fits we need further criteria in order to compare models 
based in anomalous U{1) models and non-Abelian models, such as SU{3). These 
other criteria may be found in the predictions that the models presented here can 
give in the supersymmetric sector. 
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Comparison 




U{1) (GST) 


t/(l) (Non-GST) 


5t/(3)-like 


if^ oi expansion pars. 


1 


1 


2 


^ of free pars, (quark sector) 


12 


>18 


10 


GST relation 


yes 


no 


yes 


prediction for tan (5 


small 


small 


no 


lepton sector 


o.k. 


o.k 


o.k 


simple flavour charges 


no 


yes 


yes 



Table 16: Some criteria of comparison. Here the number of free parameters corresponds 
to the number of coefHcients, phases and parameter expansions that need to be adjusted or 
determined in the fits. 

10 Flavour issues in SUSY flavour symmetry mod- 
els 

Since the flavour symmetry is expected to be broken at a high energy scale, non 
super symmetric models will have a hierarchy problem, since the cutoff of the theory 
must at least be of the order of the flavour symmetry breaking scale. Supersymmetric 
models with soft breaking parameters around the TeV scale do not have this problem. 
For this reason flavour symmetries are almost exclusively considered in the context of 
one of the minimal supersymmetric models, or one of the popular SUSY unified theo- 
ries. The soft Lagrangian parameters are strongly constrained by the supersymmetric 
fiavour problem and the supersymmetric CP problem. 

The supersymmetric flavour problem needs the soft scalar mass squared matrices 
to be diagonal to good approximation at high energy scales, since the off-diagonal 
elements contribute to one-loop flavour violating decays such as the highly constrained 
/X — > 67 in the lepton sector and 6 — > 57 in the quark sector. It also requires that 
the trilinear couplings are aligned well to the corresponding Yukawa matrix, since 
off diagonal elements in the trilinears in the mass eigenstate basis also contribute to 
highly constrained decays. The supersymmetric CP problem is related to the phases 
of the parameters in the soft Lagrangian. The general requirement is that these 
phases need to be small for the majority of soft breaking parameters. 

The reason that these problems are relevant in the context of family symmetries 
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is that in general, the existence of the family symmetry and the fields that break 
it can give dangerous contributions to the soft Lagrangian parameters. It would be 
remiss to look at these models but not check whether CP violation or fiavour violation 
is likely to rule them out. The starting point for investigating these problems is to 
consider the hidden sector part of the theory, which leads to the size and phases of 
the vevs of the fields which break the U{1) symmetry, 9 and 9. 



10.1 The flavon sector 



We start by considering the values of the expansion parameters e and e. They are 
defined by: 

MM ^ ^ 

where 9 and 9 are scalars which break the U{1)f symmetry, and have charges of 
1,-1 respectively under the symmetry. We wish to arrange that e — e, which entails 
arranging that the potential is minimized by<^>=<^>. This would be simple 
if the U{1) were non- anomalous, and thus missing a Fayet-Iliopoulis term. If we set 
the 9 sector of the superpotential to be: 

We = S{99 - Ml) (148) 

We introduce a new field, X, which has charge qx under U (1). qx will be unspecified, 
but some number such that when < X >^ 0, it doesn't contribute to the fermion 
mass operators ( or, at the very least, it doesn't contribute at leading order). Then, 
if we give 9 and 9 the same soft mass , and require that X doesn't get a soft mass, 
we end up with a hidden sector potential: 

2 

V^\99- M2|2 + ^ + _ q^\x\^ + ^-^f + rn\9'' + f). (149) 



^^This requirement may seem somewhat strong, but we also wish to minimize flavour violation 
coming from the D-term associated with C/(l), which is proportional to — m|, and will provide a 
non-universal contribution to the scalar masses. This contribution will lead to off diagonal elements 
in the SCKM basis which can easily be dangerously large with regard to flavour violation. 
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If we minimize this potential with respect to 6,6 and X, we end up with the following 
constraints: 

^ =0= 26{66-M^)+g'6{\6\^-\6\^ + qx\X\^ + e) + 2m^6 (150) 
o6 

BV _ _ _ _ 

= 0= 26{66-Ml)+g^6{\6\^-\6\^ + qx\X\^ + e) + '^rnH (151) 

o6 

_ =0= %-2X{\6\^-\6\^^qx\X\^^e) (152) 

Since X doesn't have a mass term, it would be massless unless < X >^ 0. 
Therefore, of the two solutions of Eq. ()152j) . we have to take X 7^ 0. From this, we 
see that: 

\6\^ - \W + (IxW + i'^ = ^ (153) 

Substituting Eq. p53p into Eq. p50|) and Eq. ()15H) . and multiplying by 6 and 6 
respectively, we find: 

= 66(66 - M^) + 6"^ (154) 
= 66(66 - M^) + m'^f (155) 

From this, we can deduce that either 6 = 6 = or \6\ = \6\ = Mg. The potential is 
minimized by the second solution if < 2M|. As we expect Mg to be a GUT scale 
mass, and m to be a TeV scale soft mass term, we find, that as desired, that we will 
have: 

<^>=<^>^e = e (156) 

This allows us to consider Yukawa textures without having to keep track of 
whether the overall charge for each term is positive or negative. 



10.1.1 Getting e from the Fayet-Iliopoulos term 

The GST requirement leads to needing flavon fields with opposite charges, under 
U(1)f- Were this not the case, we would have an elegant way of generating < 6 >. 
Consider a simple case where 6 doesn't have a superpotential mass term, but does 
have a soft mass: 

V = ^(-\6\^ + ef + my (157) 



52 



Then, without the need for an exphcit mass term in the superpotential, we would 
find that minimizing the potential with respect to 9 would lead to: 

<0>=i^jl + '^^i (158) 

Where the final approximation is due to the fact that we expect to be much larger 
than m?Q. So we have managed to set < > from ^, which can be predicted from 
string theory. So this allows one to predict the fiavon vev, rather than having to put 
it in by hand. 

This provides a motivation for trying to set up the case where < 6 > and < 9 > 
could both be set by the FI term. However, it doesn't seem possible to make this 
work without adding in either an extra symmetry, or extra matter. Even then, trying 
to arrange things so that < 9 >=< 9 >= zC,, with z some real number is difficult. 



10.2 Yukawa Operators 

Since the net U{1) charge can be either positive or negative and we have e = e, an 
effective potential has the following form: 

f=u,d; ij 

+ J2 L'f^Hf a^^^^+f^+'^f^. (159) 

/=e,n; ij 

We cannot say anything in particular about the Kahler potential. We can assume 
that the phases responsible for CP violation only appear in the flavour sector. Then 
observable CP violating phases will be put into the Yukawa couplings indirectly from 
the effective superpotential of Eq. p59|) . In general we can consider an effective 
Kahler potential of the form: 

K = Koit^) - ln(S + S + 5gs) + + - + Y1 ^^j"^'^' (^^0) 

i ij 

where Kq is the Kahler potential of the moduli flelds, ta = Ta + Ta, S is the dilaton, 
fi{ta) are possible functions of these moduli flelds e.g. f{t) = H^^ita^"^'^* ■ But we 
cannot specify the form of the Kahler metric. It may be that the Kahler metric is 
canonical, in which case Kf-, = 6ij*. Such a form has a good change of leading to 
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acceptable phenomenology, since the scalar mass matrices will be proportional to the 
identity at the appropriate high energy scale. When rotating the scalar mass matrices 
to the super-CKM (SCKM) basis at the high energy scale, the transformation will 
leave the mass matrices invariant. Flavour violation tends to be proportional to off- 
diagonal elements in the scalar mass matrices in the SCKM basis, so any flavour 
violation will be due to RG effects, and will therefore be suppressed. On the other 
hand, the Kahler metric could have off-diagonal structure, in which case the risk 
of flavour violating effects would be high, and the case where the Kahler metric is 
diagonal but non-universal is potentially very interesting since flavour changing effects 
are induced in general by the SCKM rotation. 

10.3 The SUSY CP problem 
10.3.1 The jjL problem 

In order to avoid the yU problem, a symmetry or other mechanism to protect /i from un- 
wanted contributions needs to be introduced. The parameter can have contributions 
from the superpotential, (expected to be at the Planck scale) and from the Kahler 
potential, via the Giudice-Masiero mechanism [11] or other mechanisms [21, 22], 
/i = nw + I^K- The charges of the flelds Hu and Hd under the flavour symmetry 
can be chosen in such a way that fiw{Mp) is forbidden in the superpotential. Then 
another fleld, S can be introduced, so that the term XSHuH^ is allowed in the Kahler 
potential, which generates an effective fi = 0{m^/2)- Note that in the cases that we 
have found for u + f 7^ there is no at Mp. In general for a theory containing 
two flavon flelds with opposite charges, once the flavour symmetry is broken below 
the Planck scale, the contributions to the term are: 

e\^+^\H^Hd^iw + e^^'+^^H^HdfiK (161) 

Thus, even if the fi term is missing from the superpotential at renormalizable level, 
it will be generated by non-renormalizable operators once the family symmetry is 
broken. However, it will appear suppressed by a factor of e'""^^'. To get an suffi- 
cient suppression, either |m + f | must be large or e must be small. Obviously, since 
the same factor el"+''| appears suppressing both superpotential and Kahler potential 
fi contributions, there is no extra constraint from considering the second term in 
Eq. (HnH). 
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However, |-u + f | is related to the anomaly cancellation conditions considered in 
Section 121 There are two possibilities for having small |m + f |. The first is to have 
small expansion parameters, e; however if e becomes too small, it makes predicting 
the fermion mass hierarchy very difficult. The second is to accept a contribution 
to n that is larger than order 0(7/13/2); however phenomenologically, the total fi 
should not be much bigger than the 0(7713/2). It is, however, possible to apply a 
new discrete symmetry to disallow the superpotential term, which never allows any 
flavon corrections to generate it. 

10.3.2 Electric dipole moment constraints 

The electric dipole moments (EDMs) constrain the form of the trilinear couplings, 
(Yj^)ij. The trilinear couplings are defined through {Yj^)ijHfQifj. Here we need to 
ensure that there is not a large contribution from the phases found in the trilinear 
terms to the CP violating phases. In the context of flavour symmetries it is usually 
postulated that the only phases appearing in the theory are in the Yukawa couplings 
and any other phase will enter as a consequence of a dependence in the Yukawa 
couplings. Then to check if the model gives contribution below the bounds one needs 
to compare the diagonal elements of the Yukawa couplings with the diagonal elements 
of the trilinear couplings, in the SCKM basis. The trilinear terms in general can be 
written as: 

(Yf^)., = Y/^F'^da [k + HKjKlK])) + F^daY/^ (162) 

We can always write the first term in a "factorisable" form [36], such that if the 
Yukawa couplings, Eq. ()159|) . are the only source of CP violation then the first term 
does not give any contribution at the leading order. For the second term, which 
involves the derivative in terms of the flavon fields, if the flavon field is the only field 
with 7^ then the diagonal trilinear couplings in the SCKM basis are real at 
leading order in the flavon fields [39]. Thus there is not an 0(1) contribution to the 
CP phases from this sector. 

One can check this simply by writing the last term of Eq. ()162j) in the SCKM 
basis: 

{F'^da{Yf)tf^^ = F^{VlMdaVLUYnr.,h + 

F^{daYBi.^)ij + F-{Yj,,,^)u{dM„{Vhr, (163) 
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Where Vl and diagonalize the Yukawa matrix: loiag = ^l^^r- The leading term 
of the Eq. ()163|) is the second term and it is at most of order 9. If another field 
has non-zero F-term, F-^ ^ then all the quantities appearing in Eq. p(j2|l can be 
written as a expansion in X and 9/M = e: 

{Yn..^)u = {au + hiX)e^^\ (164) 

We are assuming that only the matter sector in Eq. ()159p has phases leading to CP 
violation, so the term haXe^" is real and hence so is: 

F''{daY^,^^)u = FHueP^^ (165) 
10.4 SUSY flavour problem 

In addition to the F term contribution to the soft masses we have to add the D term 
contributions 

{M%, = {M^)f + {M^)d (166) 

If the Kahler metric is diagonal in the basis where the symmetry is broken both con- 
tributions are diagonal and proportional to the Kahler metric. For example, consider 
universal SUGRA: {M'^)f ij = Kijm?^. However, even if we assume that the first 
term is indeed proportional to the Kahler metric, the D-term will not in general be 
proportional to the Kahler metric: 

{M^)D^J = Y,gNXN e.K,,*{ea)ml, ml = 0{ml^) (167) 

N 

The main problems for FC processes for these kind of theories are the contributions to 
the trilinear couplings from the anomalous D-term contribution to the soft masses [26]. 
For the last issue there is no real solution so far but one can ameliorate the problem 
by making all the scalars heavier, which is a simply mass suppression. 

In order to study all the possible consequences of models with the superpotential 
structure of Eq. ()159|1 . we can parameterize the Kahler metric according to the differ- 
ent contributions it may have, assuming a broken underlying symmetry with at least 
two fiavon fields with opposite charges. Once this is done we can then study their 
consequences. As mentioned earlier, this analysis is beyond the scope of this paper. 
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so we just mention how extreme and dangerous situations may arise and we leave the 
analysis for a future reference [23] . Some authors have studied possible consequences 
of flavour models for FC effects but very specific assumptions need to be assumed 
due to the many unknown supersymmetric parameters [5,24,25]. 

The most strict bound for flavour changing processes is coming from the decay 
/i — s> 6 7 [28]- [29] and given the fact that we have a large mixing angle in the left 
handed sector of the charged lepton matrices it is crucial to determine under which 
conditions we can produce a suppressed effect. Also the constraints given by the 
process 5 — $ Ks may select out some of the possibilities presented. 

10.4.1 Non minimal sugra and diagonal Kahler metric 

Consider, for example, the case for which at the scale at which the flavour symmetry is 
broken, the Kahler metric is diagonal. For this case, we also want the soft scalar mass 
matrices diagonal but not proportional to the unit matrix, due to possible different D 
term contributions. Since the general case it is difficult to handle we consider the case 
where M| — Mf is small and M~ — M~ „ > 0. In order to estimate the flavour 
changing processes we need to take into account the effects from renormalization 
group equations (RGB's) and then at the electroweak scale make the transformation 
to the basis where the fermions are diagonal. Here we consider the case of leptons, 
since we are interested in determining 5\j and in particular which is the most 
constrained parameter due to — > e 7). 

We make an estimation of the contributions from the renormalization j3 functions 
in this case, such that at the scale where the dominant right handed neutrino it is 
decoupled we can write the soft masses as 

Mi ,(Mv.) « Ml ,(M.) - ^ in (^) ) (m 

for Mx — Mq or Mp, GUT or Planck scales respectively, and for My — Mrr 3 in this 
case and considering just one loop corrections. The (3 functions of M| from Mx 
to Mpt,R 3 receive the contributions from the MSSM particles plus the contribution 
from right-handed neutrinos. At M3 we then run from that scale to the electroweak 
symmetry breaking scale with the appropriate (3 function and matter content. In the 
case of SNRHD scenario and the form of the Yukawa matrices that we have considered 
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in Section Q we can make the following approximations for the (3 functions^^: 



xMSSM 

O -2 



L ii 



L ij 



r 2\ |2 6 2| |2 3 2o 

-6g2\m2\ - -g^\mi\ - ^Qi^ 



/ \ MSSM 

+2ml (1 + a^) (^-^^27 + ^1.3^7^7) (169) 
where we have assumed that the trilinear terms can be written as AI- = aV/.M?, 



ij e : 



and M| is not necessarily diagonal. The parameter S", defined as = — + 



Tr 



- M? - 2M? + M? + M? 



, does not generate big contributions as long the 
masses involved remain somewhat degenerate. The j3 functions generated by the 
dominant right-handed neutrino can be approximated by 



/3« J ^ 2^3^- [ml, + mUl + 6^) + m|J (170) 
From Mx = to My = Ms -the supersymmetry breaking scale-, we consider 

N MSSM 

/5^^2 j • For this estimation we ignore the effect from Ms down to the elec- 

troweak scale. At this scale we then transform the renormalized M| in the basis 
where the charged leptons are diagonal. Since there is a large mixing angle (S23) in 
the left sector of we are interested here only in estimating {Mj)^^. We can use 
the parameterization of Appendix A in order to make this transformation, i.e. 

^iag = VPY^Vi, {Mir,, = VPMlvi, (171) 

for V/j:j as parameterized in Eq. p76|) . with the j3 phases as follow 

^^For the MSSM see for example [27], when including right handed neutrinos, see for example [28]. 
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Using these approximations, we obtain the following results 



{Ml)[2 = Sl2(c2>i 22-^^111) + 

K)'i3 = c^^.^t^^§e^"3L(^| ^^_^2 

+ t23e*^m| ^2 + "^1 13) . 
("^1)23 = ^235236*^^^ (m| 22 — ^\ 33) 

+e^"-cS ((c^2t)'e-^^m| 23 - (^2t)'e^^/5;. ^ j , 

K)n = («l§)'((c2t)'^i22 + (^2t)K33) 

K)22 = {c^f {{c^^fK 22 + (^2§)M 33) - (cS)'c-.-2Re{m? 236"^^ 

K)33 = (cS)M(^2t)M22 + (c2t)'^i33)+C2t^2n2Re{m?23e-'n) 

(173) 

here the soft masses are the soft masses at Ms, renormalized from Mx = Mq, Mp 
down to M3 with the appropriate contributions from the dominant right handed 
neutrino, Eq. (IT!)Hll . and Eq. ^T^ -Ea. (IT7n|l and then from M3 to Mg with the 
appropriate P^MSSM) functions. Thus we began with a diagonal matrix M| at Mx, 
then the RGE effects up to the scale where M3 is decoupled generate a non diagonal 
matrix which receives more RGE contributions from M3 to M5. At electroweak scale 
we transformed to the basis where charged leptons are diagonal. The mixing angles 
in this sector can be approximated as 

42 = l(«12 - ^32013) I/I («22 - «32«23)|e''^^ ^^3 = a?3/«33e''"> ^3 = «23/«33 (174) 

The powers pf^- for the different solutions presented now correspond to = 2/3, 14/3, 
pI^ = 29/12, 71/12 for Fits 2 and 3 respectively. So in this case we see that we need a 
big suppression of the element (m? in order to be in agreement with the observed 
bound on /i — > 67. In the present example the suppression it is related to a bound 
on ijn\i~f^\2) s^nd a relative big set of soft masses. The results of these estimations 
are presented in Table IT7I 

As we can see from the results of Table (fT7|) the estimation of |(5ii)i^| is less 
dependent on the relation among the original soft mass terms than on the value 
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Estimation of 6ij for the Fit. 3 of Section IHl 


Paramter 


Ex. I 


Ex. II 


Ex. Ill 


m^JGeV] 

^L3[GeV] 
rriei [Ge V] 
me2[GeV] 
mesiGeVJ 
Ml [GeVJ 
M2[GeV] 
MnAGeV] 
M^JGeV] 
Ms 


520 
530 
500 
520 
530 
500 
500 
2Mi 
510 
510 
1000 


520 
530 
500 
520 
530 
300 
500 
700 
510 
510 
1000 


520 
570 
230 
520 
550 
300 
500 
700 
510 
510 
1000 




514 


486 


456 


2 / 2 

X = mi^/mj 

I&)f2l 
I&)f2l 


4.3 X 10^3 
0(1 
1.7 X 10-3 
5.7 X 10-2 
0(1 


0.3 

5.6 X 10-3 

0-1) 

1.8 X 10-3 
6.4 X 10-2 

0-1) 


1.4 X 10-3 
0(10-2) 
1.9 X 10-2 
6.3 X 10-1 

o(io-i) 



Table 17: Estimation of in the fit 3 presented for the non minimal sugra example 
and its comparison to the observed bounds \6lj\^ [28]- [29]. 

taken for the average s-lepton mass, which indeed needs to be large. Here we note 
that this is just an estimation on the conditions that B{fi —^ej) imposes on the 
soft masses, but with out fully checking whether or not appropriate masses for all 
the MSSM parameters can be obtained. In the following we consider a numerical 
investigation in the minimal sugra case. 

10.4.2 Numerical Investigation of B{fi ^ e 7) in minimal sugra 

The presence of a right-handed neutrino fields leads to RG lepton flavour violation. 
Since the masses of the right handed neutrinos are so light for the GST solutions, fits 
1-3, we attempted a numerical analysis for all of the fits of Section (0) using the same 
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modified version of SOFTSUSY [41] as used in [37]. 

In order to get a good liandle, we liave embedded tlie flavour model fits into a 
string-inspired mSUGRA type scenario, with no D-term contribution to the scalar 
masses. This scenario was chosen because it is expected to be the embedding with 



mass 7713/2. 

As ni was only constrained to be between —ajl and 0, we allow it to vary within 
this range. We define the model at the GUT scale as: 



This setup of the soft parameters corresponds to benchmark point A in [37]. The 
results are as follows, for Fit 1 the code being used can not generate any low energy 
data for this fit so we do not find any safe B{fi —>■ ej) region using the conditions 
presented above. The Fit 2 has BR(/i — > 67) <= 10~^° which is unattainably low, 
thus this fit is plausible within the context of the minimal sugra conditions that have 
been specified. The smallness of the branching ratio for fit 2 comes about because 
with no RG running, in mSUGRA this rate would be exactly zero. The RG flavour 
violation will come from terms proportional to V^^V^, whose elements are tiny ( the 
largest is 0(10^^^) ). 

The Fit 3 generates a tachyonic s-electron for the full (7713/2,^1) range. This is 
not to say that this fit will always have a tachyonic s-electron in other, less trivial 
embeddings. Fits 4 and 5 produce regions below and above the experimental limits 
on B{fi —>■ 67), the graphs for these fits appear in Tables (00)). 

11 Conclusions 

In summary, we began our analysis by reviewing the Green- Schwartz (GS) conditions 
for anomaly cancellation for theories based on a U{1) family symmetry. We then 
used these conditions to fix the charges of all the quark, lepton and Higgs fields and 
studied possibilities where the Higgs mass /i term is either present or absent in the 
original superpotential. The solutions which we constructed do not necessarily require 
an underlying Grand Unified Theory (GUT) but may be consistent with unification 
because of the GS conditions. Regardless of the presence of an explicit unified gauge 



the lowest flavour violation. In the scenario, Aq, 



,Q, Ml/2 are all related to a gravitino 



m. 




(175) 
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Figure 1: BR(/i ^ 67) for fit 4, with (S) = 0(Mg). The sohd points are below the 
experimental limit of 1.1 ■ 10^^^, and the hollow points are above. 
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Figure 2: BR(/i — 67) for fit 4, with (S) = 0{Mpi). The solid points are below the 
experimental limit of 1.1 ■ 10^^^, and the hollow points are above 
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Figure 3: BR(/i ^ 67) for fit 5, with (S) = 0(Mg). The sohd points are below the 
experimental limit of 1.1 ■ 10^^^, and the hollow points are above 
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Figure 4: BR(/i — 67) for fit 5, with (S) = 0(Mp;). The solid points are below the 
experimental limit of 1.1 ■ 10^^^, and the hollow points are above 
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group, the explicit solutions can produce matrices of the form that are identical to 
those that would be expected in an SU{5) case or Pati-Salam unified theory, for 
example. 

The flavour structure of the resulting Yukawa matrices is controlled by the charges 
of the quarks and leptons under the U{1) family symmetry gauge group. We have 
determined these charges which are consistent with anomaly cancellation, and studied 
cases which can reproduce quark Yukawa matrices satisfying the Gatto-Sartori-Tonin 
(GST) relation, as well as other cases which do not satisfy the GST relation. We 
find the GST relation to be an appealing description of the value of the element Vus, 
and the GST relation provides a useful criterion for classifying flavour models. In our 
view, having the Cabibbo angle emerging automatically from a flavour model should 
have a similar status to gauge coupling unification in a high scale model. Having 
classified the solutions in terms of the GST condition, we then further classify the 
solutions according to which of them can produce the observed mixings in the lepton 
sector, and those that are consistent with a sub-class of solutions based on the SRHND 
or sequential dominance scenario with the further condition that the charges of the 
lepton doublets for the second and third family are equal, h — h- We find that the 
GST solutions combined with SRHND results in highly fractional charges. On the 
other hand non-GST solutions with SRHND results in simpler charges, and we have 
therefore studied both sorts of examples. 

We have presented three numerical examples of solutions satisfying the GST re- 
lation and two examples of non-GST solutions in order to compare how well these 
solutions fit the experimental information while maintaining 0(1) coefficients. For 
the GST solutions, one of these examples corresponds to a model that can be thought 
of as coming from an underlying SU (5) and for which a ^l term is allowed in the su- 
perpotential. It is well known that in this case, given the relation — ^, there 
should be a Clebsch-Gordan coefficient different in the charged lepton (2, 3) sector 
and in the (2, 3) d-quark sector in order to produce appropriate mixings in the con- 
text of the C/(l) flavour symmetry and the GUT theory. Two other GST examples 
are presented for which the term is not allowed and which are not consistent with 
an underlying SU{5), or other GUT theory. In these cases ^ {Y^Y but it is 
possible to maintain the relation nib and in one of them just the 0{1) coeffi- 

cients of the underlying U{1) theory can account for the appropriate mixings in the 
charged lepton and d-quark sector. The non-GST cases also give a good description 
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of masses and mixings, although in this case we need to rely on further coefficients, 
possible Clebsch-Gordan coefficients from an underlying GUT, in order to achieve a 
good phenomenological description. 

For the above examples we have provided detailed numerical fits of the 0(1) 
coefficients required to reproduce the observed masses and mixings in both quark 
and lepton sectors. The purpose of performing such fits is to compare how well the 
different models can fit the data, and to try to determine quantitatively the best 
possible model corresponding to the best possible fit. Although in the cases just 
mentioned the solutions which fit the data best are the solutions consistent with 
an underlying SU (5) theory, the other two fits are quite plausible and represent 
interesting possibilities which cannot be excluded. Since all the models constructed 
have good agreement with the fermion masses and mixings, we clearly need further 
criteria in order to discriminate between the different classes oiU{V) family symmetry 
models. 

One may ask the more general question whether family symmetries based on 
abelian or non-abelian gauge groups are generically preferred? In order to address 
this question, we have extended the fit to include a generic symmetric form of quark 
and lepton mass matrices that can be understood in the context of a theory based 
on SU (3) family symmetry. We have found that overall the generic SU (3) family 
symmetry produces Yukawa matrices which tend to fit the data better, although the 
effect is not decisive, and one cannot draw a strong conclusion based solely on fits to 
fermion masses and mixings (or the way they can be reproduced). We have there- 
fore enumerated some other possible criteria that are important in order to further 
discriminate among different flavour theories. Including the effects from the super- 
symmetric sector provides an additional way to discriminate among different theories 
based on their different predictions for soft masses and the resulting flavour changing 
processes and CP violation. We have presented two frameworks in which these pro- 
cesses can be studied in the context of flavour theories. The flrst is a non-minimal 
sugra scenario where family symmetries may render the Kahler metric diagonal at 
the flavour symmetry breaking scale, with off-diagonal elements arising only due to 
RG contributions and the non-degeneracy of soft masses. The second framework is a 
minimal sugra scenario for which a numerical exploration of — > e 7 was performed. 
The results of this analysis shows marked differences between the different models 
presented. Of the GST cases only one survives the test of S(// — > e 7) while for all 
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of the non-GST cases presented there exist regions compatible with the B{ij, ^ e 
experimental limit. 

In conclusion, at the present time, phenomenological analyses provide some guid- 
ance about what family symmetry approaches may be valid, but do not yet allow one 
to draw any firm conclusion. More specific assumptions or data in the supersymmet- 
ric sector arc needed in order to further discriminate between classes of models based 
on different family symmetry, unification or GST criteria. 
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A Conventions for the Yukawa diagonalization ma- 
trices 



We diagonalize the Yukawa matrices, , with the unitary matrices and such 
that l^Diag = VpY^V^^. can be parameterized as 



(176) 
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where a plane R^ rotation has the form: 
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°2Z 


•-23 



(177) 



In this notation, the CKM matrix is F = K'V/ 



B Comparison to experimental information 

The experimental information determining Vckm, usually put in terms of the Wolfen- 
stein parameters A, A, p and r], is extracted mainly from semileptonic decays of B 
mesons, CP violation in the K system, B'^^ — B^^ oscillations, and CP asymmetries 
in various B decays. We use a fit, based in a Bayesian approach (see for example [19] 
and [40]), of the the parameters A, A, p and rj including all the available information. 
Once we have done this, we compare the predictions of the mass textures with the fit- 
ted parameters because these include in a statistical way the experimental information 
from all the experiments considered. In the limit where we neglect all super symmetric 
contributions to these observables, the fitted values for p and fj are 

p = 0.199t°:°i, fj = 0.328l°:°i6^ (178) 

In order to compare to the Vckm prediction, as given by the U{1) symmetries, we 
need to choose four elements, or combinations of them, we choose 

wy m '"f^*' ^"'^ 
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where J is the Jarlskog invariant. We choose these parameters because they can be 
put neatly in terms of the Wolfenstein parameters 



IK 



ub 



cb\ 



A 



Wt 



td\ 



\Vts\ 

Im{J} 



IK. 



A, 



To include the information of the quark masses we use 



fUr. 



rrir. 



mt 



rrih 



:i80) 



;i8i) 



The ratios — and — can be determined from the best measured ratios of the fol- 
lowing mass ratios and the Q parameter, which is determined accurately from chiral 
perturbation theory; 

nxs/nid 



rrid 



trie 
rris 



Q 



'rriu \2 



(182) 



We note here that a change in rris with respect to previous similar fits [19] has an 
impact in the coefficients determined for the SU{3) symmetry, although it is consistent 
with previous determinations if we consider the errors involved. We have used here 
mc/rris = 15.5 ± 3.7 in contrast to mc/rris = 9.5 ± 1.7 as used in [19]. We put in 



Parameter 


Exp. value 


Value at Mx 




(9.16 ±0.67) X 10-2 




Vtd 

Vts 


0.1989 ±0.0093 




Vus 


0.224 ±0.0036 




Im{J} 


(2.88 ±0.4) X 10-5 


(1.4 ±0.5) X 10-5 


ruu 
nic 


(1.9 ±0.19) X 10-3 




rUc 
mt 


(7.5 ± 1.7) X 10-3 


(2.6 ± 1.8) X 10-3 


"M 
rus 


(5.2 ±0.35) X 10-2 




rris 

■nih 


(6.4 ±2.3) X 10-2 


(4.5 ±2.4) X 10-2 



Table 18: Experimental values used for the fit of numerical coefficients in the quark 
Yukawa matrices 

Table ()18|1 the experimental values of the parameters that we use to determine the 
coefficients of the Yukawa texture. From equations Eq. (jl79|) and Eq. (jl8ip we see 
that we can fit only eight parameters of both of the Yukawa matrices as given by 
the U{1) symmetries, but we will see that in most cases that is sufficient in order 
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to account for the viability of a given ansatz or symmetry for the matrices. We also 
need to take into account the RGE effects in going from the electroweak scale to the 
scale at which the U{1) symmetry breaks. We assume first that this is the GUT scale, 
so in order to determine the values of the parameters defining the U{1) symmetry 
(analogously for the SU (3) case) we take the values of the parameters appearing in 
Eqs. fjl79H18l|) at the GUT scale. One of the reasons for using mass ratios instead of 
just masses is because the RGE effects on the mass ratios has less impact than for the 
masses. This fit of the parameters defining the U{1) symmetry is performed with the 
aid of the MINUIT package adapted for root [42] . In this way we are able to compare 
how well a symmetry is fitted to the experimental information, and compare among 
the fits for different symmetries. 

We do a completely analogous analysis in the neutrino sector, using the following 
observables 

in 

tan 6*23, tan ^13, tan ^12, m,^,, — — (183) 
and their experimental values as appear in Table (|T9|) 



Parameter 


Exp. value 


tan ^23 


1.07 ±0.37 


tan^^i3 


0.21 ±0.1 (u.b) 


tan 6*12 


0.65 ±0.12 




0.05 ±0.01 


"1^3 


0.19 ±0.05 



Table 19: Experimental values [43] used for the fit of numerical coefficients in the 
neutrino Yukawa matrix. For tan ^13 we have fitted using the upper bound. 

B.l Evaluation of observables with fitted parameters 

In this section we put the evaluation, at the scale Mx, of the experimental inputs 
using the fitted parameters in order to compare with Table ^1 
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GST sol. 2u =t; = 


GST sol. 2u^v^ 


GST sol. 3 u^v^O 


Parameter 


Evaluation at Mx 


V h 

' CO 


(9.15 ±0.52) X 10-2 


(9.12 ±0.35) X 10-2 


(9.31 ±0.27) X 10-2 


yts 


0.199 ±0.081 


0.199 ±0.031 


0.20 ±0.078 




0.225 ±0.0041 


0.225 ±0.0043 


0.225 ± 0.0040 


Im{J} 


(1.4 ±0.3) X 10-5 


(1.4 ±0.3) X 10-5 


(0.8 ±0.6) X 10-5 


TTlc 


(1.9 ±0.21) X 10-3 


(1.9 ±0.16) X 10-3 


(1.8 ±0.21) X 10-3 


rric 


(2.1 ±0.73) X 10-3 


(1.9 ±0.47) X 10-3 


(4. ± 1.54) X 10-3 


mg 
nis 


(4.9 ±0.5) X 10-2 


(4.9 ±0.5) X 10-2 


(4.9 ±0.6) X 10-2 


ms 
rrih 


(3.0 ± 1.7) X 10-2 


(1.9 ±2.1) X 10-2 


(4.4 ±1.3) X 10-2 


tan 6^23 


1.27 ±0.39 


1.02 ±0.41 


1.00 ±0.42 


tan ^13 


0.21 ±0.32 


0.21 ±0.39 


0.27±0.19 


tan^^i2 


0.65 ±0.41 


0.65 ± 0.42 


0.65 ±0.32 




0.046 ±0.05 


0.023 ±0.05 


0.02 ±0.05 


—A 


0.19 ±0.05 


0.16 ±0.08 


0.19 ±0.05 




Non GST sol. lu=v=Q 


Non GST sol. 2u =v =0 




Parameter 


Evaluation at Mx 


V h 

^ CO 


(9.59 ±0.49) X 10-2 


(9.14 ±0.65) X 10-2 




Vtd 

yts 


0.201 ±0.073 


0.199 ±0.033 




Vus 


0.225 ±0.0037 


0.225 ±0.0027 




Im{J} 


(1.7 ±0.4) X 10-5 


(1.4 ±0.3) X 10-5 




rriu 

TTlc 


(1.9 ±0.33) X 10-3 


(1.8 ±0.15) X 10-3 




roc 

TtXt 


(1.1 ±0.82) X 10-3 


(2.5 ±0.56) X 10-3 




nis 


(5.2 ±0.3) X 10-2 


(4.9 ±0.5) X 10-2 




rris 
rrih 


(3.2 ± 1.2) X 10-2 


(4.5 ±2.2) X 10-2 




tan 6123 


0.67 ±0.41 






tan 6*13 


0.19 ±0.16 






tan 6^12 


0.55 ±0.31 








0.038 ±0.04 






m^3 


0.13 ±0.06 







Table 20: Evaluated values of experimental observables using the fitted parameters 
for the SU (5) GST case, for the u ^ —v ^ Q cases, the SU{b) non GST case and the 
S'f/(3) like cases considered in Section 
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SU{Z), $1 < 


SU{Z), $1 > 


Parameter 


Evaluation at Mx 




(9.14 ±0.65) X 10-2 


9.13 ±0.21) X 10-2 


Vtd 

Vts 


0.199 ±0.023 


0.198 ±0.024 




0.225 ± 0.0027 


0.225 ± 0.0033 


Im{J} 


(1.4 ±0.3) X 10-5 


(1.4 ±0.3) X 10-5 


rriu 
trie 


(1.8 ±0.15) X 10-3 


(1.9 ±0.15) X 10-3 


nic 
nit 


(2.5 ± 0.56) X 10-3 


(2.5 ±0.56) X 10-3 


rrid 
rris 


(4.9 ± 0.5) X 10-2 


(4.9 ± 0.5) X 10-2 


rris 
nih 


(4.5 ± 2.2) X 10-2 


(4.5 ±2.2) X 10-2 



Table 21: Evaluated values of experimental observables using the fitted parameters 
for the SU (3) like cases considered. 
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